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A pure sample of CH;D has been prepared and studied in the mass spectrometer. The mass 
spectrum of CH;D cannot be predicted from that of CH, by assuming complete similarity 
between H and D (except mass). More information has been obtained on the difference in 
behavior of H and D on electron impact previously noticed by Delfosse and Hipple for CzH2D>2. 
There is a difference in dissociation probability of C—H and C—D bonds which is approxi- 
mately constant for several ionization processes and over a considerable range in electron 
energy. In addition to this there is a difference in behavior of H atoms in CH, and H atoms 
in CH;D. In order to account for the mass spectrum of CH;D in terms of the mass spectrum of 
CH, both these effects must be taken into consideration. Re-examination of the data of 
Delfosse and Hipple points not only to the difference between H and D atoms in C2H2Ds2, but 
also to a difference in behavior of H atoms in C,H, and H atoms in C2H2Ds». 





INTRODUCTION 


ELFOSSE and Hipple,' in a mass spectro- 
metric study of di-deutero-ethylene, cis and 
trans, showed that upon dissociation by electron 
impact the splitting of the C—H bond is favored 
over the splitting of the C—D bond. In ac- 
counting for ions of high mass, which are formed 
by breaking C—H and C—D bonds only, they 
obtained a factor of 0.3 for the relative proba- 
bility of dissociation of a C—D and C—H bond. 
For ions in the low mass region, formed by 
breaking the carbon-carbon bond in addition to 
carbon-hydrogen bonds, a factor of 0.5 was ob- 
tained. They also compared the relative abun- 
dance of corresponding ions in C2H, and C,H2Dz, 
and found that they were nearly the same with 
the exception of the (C:H;+) type. We believe 
this exception to be significant. 


1J. Delfosse and J. A. Hipple, Phys. Rev. 54, 1060 
(1938). 


In an effort to reduce the complicating factors 
of interpretation to a minimum we have prepared 
the simpler molecule mono-deutero-methane and 
obtained its mass spectrum. This spectrum has 
been examined with a view to determining if the 
selective process observed for di-deutero-ethylene 
is a general phenomenon. Secondly, the spectrum 
for CH3D has been compared to that for CH, 
in an attempt to discover the reason for the 
anomalous behavior of the (C2:H;+) type ion 
from CoHeDo. 

We obtained spectra for a series of electron 
energies ranging from 55 volts down to the region 
of the appearance potentials of the several ions, 
At the lower energies the effect of ions formed by 
the liberation of more than one hydrogen type 
atom is reduced to a minimum. Furthermore, 
the generality of the phenomenon of selective 
splitting is tested over a wide range of impact 
energies. 
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TABLE I. Relative intensities of ions of given m/e for 
CH;D and CH, produced by 55-volt electrons. 











Ions CH;:D CH, 
(CH,*) 100.0 
(CH;D*) 100.0 
(CH2D*) 75.6 +0.6 
(CH;++CHD?) 22.4 +0.1 82.6 
(CH.*++CD*) 9.41+0.10 17.0 
(CH*) 5.20-+0.07 8.19 

+ 2.41+0.03 2.35 
(HD*) 0.12 
(H2++D*) 0.74+0.02 0.35 
(H+) 4.60+0.10 5.45 








Finally, it is of interest to test the existence 
of the effect for the splitting off of Ht and D* ions 
as well as of H and D neutral atoms. 


EXPERIMENTAL 


The purity of the CH;D, particularly with 
regard to the presence of CH,4, was of utmost 
importance. The compound was prepared by the 
reaction of D,O with the Grignard reagent 
methyl magnesium iodide in diethyl ether. The 
total volume of the reaction unit was 75 cc. 
Hence, the relative amount of H,O remaining on 
the glass after forty-eight hours of heating at 
100°C should be negligible in comparison to the 
one cc of D.O used. In order to test possible 
water contamination during assembly of appa- 
ratus and transfer of reagents the vapor in 
equilibrium with the Grignard solution was 
withdrawn at reduced pressure and analyzed 
for CH, Within the limit of detection (0.1 
percent) no CH, was found. 

The product CH;D was fractionated from the 
ether until no ether peaks were observed in the 
mass spectrum. The final product showed no 
m/e peak greater than 17. 

The D,.O was analyzed both by density meas- 
urements and by conversion to D2 with subse- 
quent analysis in the mass spectrometer. Both 
analyses showed that the sample was 99 percent 
D.O or better. 

The mass spectrometer was the Consolidated 
Engineering Corporation instrument.? The gas 
pressures used were low enough to insure that 
only the primary products of electron impact 
were observed. 


2H. W. Washburn, H. F. Wiley, and S. M. Rock, Ind. 
Eng. Chem. Anal. Ed. 15, 541 (1943). 


For the ionizing efficiency curves the nominal 
values of electron energy, as set on the instru- 
ment voltmeter to +0.3 volt, were corrected to 
include the effect of ion draw-out or pusher 
voltages. Although for each different m/e the 
ionizing electrons acquire different amounts of 
energy from the pusher voltage because the 
pusher voltage is proportionate to the ion ac- 
celerating voltage, the proper correction at each 
m/e can be derived easily from calibrations with 
a gas of known appearance potential. By calibra- 
tions with argon, the ratio between electron 
energy correction and ion accelerating voltage 
was determined with an accuracy of about 5 per- 


_ cent, allowing the assignment of electron energies 


to within 0.5 volt. The difference between the 
appearance potential of CH;D ions and that of 
CH, ions is influenced only slightly (~0.1 volt) 
by the probable error in the pusher voltage 
corrections. Decrease of gas pressure during the 
measurement of ion intensity as a function of 
electron energy was minimized by the use of a 
large sample reservoir (15 liters), and a small 
correction for this decrease was made. 

Three samples of CH;D were prepared, and 
satisfactory agreement of the spectra was ob- 
tained. For purposes of comparison a sample of 
CH, was prepared in the same manner. 


RESULTS 


The average values obtained for the relative 
ion intensities for CH;D are listed in Table I, 
together with those obtained for CH, under 
similar instrument conditions. The ion intensity 
per unit pressure is the same for CH;D* and 
CH,"*, within 1 percent. 

A simple relationship might be expected be- 
tween the mass spectra of CH3;D and CHy,. One 
can try to predict the spectrum of CH;D by 


TABLE I]. Spectrum of CH;D predicted from that of CH, 
assuming identical behavior of H and D. 











CHs Predict 

Ion spectrum CH;:D 
CH;D* 100.0 

CH,* 100.0 

CH:Dt 62.0 
Pathe 82.6 20.6 
D 8.5 
CH,* 17.0 8.5 

CDtr 2.04 

CH+t 8.19 6.15 
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assuming (1) that H and D atoms behave in 
identical fashion in CH;D, and (2) that H’s in 
CH;D and CH, are likewise identical. It is 
furthermore reasonable to assign ion intensities 
for the CH;D spectrum by applying simple 
statistics to the CH, spectrum; e.g. (CH;* from 
CH;D) =3(CH;* from CH,). In this way the 
values given in Table II are obtained. In Table 
III the values of column three, Table II, are 
summed in order to compare the predicted in- 
tensities of ions for a given m/e with the intensi- 
ties actually observed. The disagreement in 
Table III shows that the above simple approach 
is not adequate. 

We proceed to test our first assumption, 
namely, that H and D are identical in deutero- 
methane. With this assumption we shall find 
that the apportioning of intensities observed at 
each m/e to individual ions of the CHs;D spec- 
trum results in internal inconsistencies, as shown 
below. Equations (1)—(3) express the relation- 
ships to be tested. 


(CH2D*) =3(CH3*), (1) 
3(CD+*) =(CH?*), (2) 
(CHD*) =(CH2*). (3) 


Substituting the appropriate values of CH;D 
peaks (Table I), we have: 


(CHD*) = 22.4—(CH3+) =22.4—75.6/3 = —2.8 
and 
(CH2*) =9.41 —(CD*) = 9.41 —5.20/3 = 7.68. 


Consequently (CHD*) #(CH;*) and assumption 
(1) is therefore incorrect. We conclude that the 
same selectivity of splitting shown with C,H.D» 
is exhibited here. 

It should be possible to calculate the relative 
probability of breaking a C—D bond, and the 


TABLE III. Observed spectrum of CH;D compared with 
predicted spectrum assuming identical behavior of H 
and D. 











m/e Observed Predicted from CHs 
16 75.6 62.0 

15 22.4 29.1 

14 9.41 10.54 

13 5.20 6.15 








TABLE IV. Relative intensities for individual ions of CH;D; 
predictions from CH, using rp =0.38. 











I SD Cin 
on 3 3 
CH;D+* 100.0 100.0 
CH:D+ 75.6 62.0 
CH;* 9.6 7.8 
CHD+ 12:8 8.5 
CH;+ 4.85 3.2 
Cp* 4.56 2.05 
CH* 5.20 2.34 








correlative factor can be defined as follows: 
probability of D leaving CH;D 
~ probability of individual H leaving CH;D’ 





TD 


We make the simplifying assumption ‘that 7p is 
constant for the formation of all types of ions 
from CH;D. Equations (1)—(3) are now written 


3(CH3*) = mp(CH2D*) 
or (4) 
(CHD*) = 22.4—75.62p/3, 


(CH+) = 3ap(CD*) 
or (5) 
(CH;+) =9.41—5.20/3mp, 


(CH2*) =mp(CHD*). (6) 


Solving, we find that mp for deutero-methane 
when 55-volt electrons are used is 0.38, a figure 
with which the Delfosse and Hipple value of 
0.3-0.5 compares reasonably well. 

The effect of subtraction from the ion intensi- 
ties of the contribution of 1 percent CH, due to 
possible H,O in the D.O is to lower mp to 0.37, 
which is within the experimental error. 

We can also make an estimate of the numerical 
value of the mp* factor, which will express the 
relative ease of breaking C—D and C—H bonds 
to form D* and H* ions. 

We can place an upper limit on this value by 
assuming that no H,* is formed in the ionization 
chamber, thus attributing the entire 2 peak 
(0.74) to the D* ion, and finding 


3(0.74) ; 
48. 


1p*t (max) =————- = 0. 


The true value will be smaller than 0.48, for 
some H,* is undoubtedly formed, and this would 
make the intensity attributable to D* smaller. 
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TABLE V. “Observed” and predicted CH3D spectra assum- 
ing 12 percent methane impurity. 


BAUER, AND BEACH 


TABLE VI. Comparison of intensities of corresponding ions 
of CH;D and CH,. 

















; ae Predict y CH, 
on serve Tp =0. 
CH:Dt+ (62.0) (62.0) 
CH;* 5.18 5.15 
CHD+ 6.02 8.5 
CH,* bee py 
CDt 5.48 2.05 
CHt 4.12 1.54 








The CH, pattern, for example, shows an (H2*) 
peak of 0.35. If one-half as much (H2*) is formed 
from CH3;3D as from CH,, 


3 
mpt =——(0.74—0.5 X0.35) =0.37. 
4.60 


A check on this second value is to assume that 
the sum of the corresponding ions H2+ and HD* 
in CH3D equals H2* in CH,y. Then 


and 


3 
rpt =——(0.74—0.23) =0.33. 
4.60 


The best that can be said is that zp* is certainly 
less than 0.5 and probably greater than 0.3, an 
interesting result in view of the value mp =0.38. 

We can now return to the problem of finding 
a correlation between the CH;D and CH, 
spectra, correcting our first assumption by intro- 
ducing 7p in the following way. With the known 
value of zp applied to Eqs. (4)—(6), the observed 
peak heights of CH;D (Table 1) can be broken 
down into ‘‘observed”’ contributions due to indi- 
vidual ions, shown in column two of Table IV. 
Similarly, the value of mp is used to predict 
(Table IV, column three) the individual ion 
intensities of CH;D from the CH, spectrum; 
e.g.,-(CH;+ from -CH;D)=2p/4(CH;* from 
CH,). 

The predictions and observations are still at 
wide variance, though all differences are now in 
the same direction. The first thought is to 
attribute the disagreement to a CH, impurity, 
despite precautions taken against this occurrence. 
This explanation cannot suffice even by attrib- 
uting to methane all of the discrepancy observed 
for the 16 peak, which amounts to 75.6—62.0 





Ion type CH;:D CHa 
CX,t 100.0 100.0 
CX;t 85.2 82.6 
CX_t 17.65 17.0 
CxXt 9.76 8.19 
co 2.41 2.59 








= 13.6 units. This amount of CH, (12 percent) 
would necessitate a corresponding correction to 
the CH;D spectrum and would cause mp to be 
lowered from 0.38 to 0.25. Table V compares 
the predicted and ‘‘observed’’ spectra on the 
basis of tp =0.25, and shows that this extreme 
correction fails to remove all the discrepancies. 
Furthermore, the fact that the ion intensity per 
unit pressure was found to be equal within 1 
percent for CH,* and CH;D* supports the ex- 
pectation that both the CH, and the CH;D 
preparations contained no more than about 1 
percent impurity. It is reasonable, therefore, to 
regard the differences indicated in Table IV as 
real and caused by the incorrectness of our 
second original assumption, namely, that H’s in 
CH;D and CH, are identical. In fact, we shall 
see that these differences mean that the presence 
of the D atom in the molecule causes the H 
atom to be removed more easily from CH;D 
than from CH,. 
We define a new quantity In. 


probability of individual H leaving CH;D 
Tqy= 





probability of individual H leaving CH, 


TABLE VII. Relative ion intensities for C2H, and C2H2D2. 











Predict Predict ry 
Observed Observed C2H2D2 CeH2De 
Ion C2H2D2 CoH Tp =1 zrp=0.33 mp =0.33 

C:H2D.2+ 100 100 100 100 
C:H.Dt 11 60 30 10 

‘33 
C,:HD,* 34 30 30 
C.D.* 24 9.5 9.5 
C.HDt 32.5 38 12.6 1.59 
C.H.* 2.6 57 9.5 - 1.03 
C.D+t 8.2 4.8 1.59 

1.73 
C:H* We | 9.6 4.8 0.53 








® See reference 1 in text. 





Th 
CE 


ref 


in | 


r'x( 


For 


El 
(v 





gz ions 


shall 
ence 
ie H 
H;D 


H;D 
‘H, . 





H2Dsz. 








, =0.33 





1.59 


MASS SPECTROMETRIC STUDY OF CH;D 705 





+ 
CH FROM CH, 










CH,D* FROM CH4D 


CH,* FROM CH,D 


oe —_ 4 





| 1 i j — 





90r- 

80F 

ta 7Or 
= 

7 60b 
Fic. 1. Voltage of the ionizing elec- 2 

trons versus relative ion intensities of © sol 
CX;+ type ions. The CH;* ion from z 
methane is relative to the CH,*+ ion; ~ 

CH;*+ and CH:D* ions from mono- z “9 
deutero-methane are relative to the © 

CH;D* ion. Ionizing current is con- 30F 
stant. 4 
> 

p 20 
< 

GJ io 
fo 

= 
Oo i2 14 


The subscript DM will refer to the compound 
CH;D, M to CHy. The expression (CH2D+)x 
refers to a hypothetical ion formed from CH;D 
in which H and D both behaye exactly like H 
in CH,. 


(CH2D*)m =3/4(CHs*)m, 
T'n(CH2D*)m = (CH2D*+)pm, 
(CH2D*)pm 

(CH3*)s 
=4/3X75.6/82.6=1.22 


Ig = 4/3 x 


For the ions (CH2*)pm and (CHD*)pm 
l'g?(CHD*)m = (CHD*)pn, 
ashes 2(CHD*)pm < 2X12.8 
(CHe*)m 17.0 
ly = 1.23. 





= 1.507, 


TABLE VIII. rp for CH3D over a range of 
electron energies. 














Electron Electron 
energy energy 
(volts) Tp (volts) Tp 

55 0.38 23 0.41 
32 0.40 22 0.40 
30 0.40 21 0.42 
29 0.40 20 0.42 
28 0.41 19 0.43 
27 0.40 18 0.44 
26 0.40 17 0.44-0.46 
so 0.39 16 0.35-0.47 
24 0.40 15 0.69-0.43 








= —— 
22 22 24 26 28 30 32 34 36 38 


IONIZING VOLTAGE 


For the ions (CH*+)pm and (CD*)pm 
l'p?(CD*)m - (CD*)pm, 


4X4.56\! 
ra~(——) ==] 31 
8.19 


It is apparent that two effects exist. The first 
is the “x effect,’’ which measures the difference 
in activity of an H and a D on the same carbon 
atom. This amounts to a factor of about 2.6. 
The second is the “‘T effect,’’ which is a specific 
activation of H atoms on CH;D over H on CH,. 
The hydrogen atoms on deutero-methane are 
removed 20 percent—30 percent more easily than 
those on ordinary methane. 

The fact that the sums of the CH;D corre- 
sponding ions agree fairly well with the corre- 
sponding ions of CH, (Table VI) appears to be 
fortuitous and the result of two opposing effects 
which in large part cancel each other. 

Delfosse and Hipple found similar rough 
agreement of intensities of corresponding ions of 


’ CsH2D2 and C2H, with the exception of the ion 


type (C2X;+). In this connection we have re- 
examined their data to see if this discrepancy 
can be removed by including a gamma effect. 
Table VII shows the C2H,and C2H2Dz spectra, 
with the spectrum for C:H2Dz predicted from 
C.H, using rp =1 and rp =0.33. The last column 
gives the values of 'y which we have calculated 
for different ion types of the molecule. The 
values do not show the same constancy observed 
in the case of CH;D, but it is apparent that there 














CHD* FROM CH3D 


CH3 FROM CH,D 


BAUER, AND BEACH 


CH} FROM CH, 


Fic. 2. Voltage of the ionizing elec- 
trons versus relative ion intensities of 
CX.* type ions. The CH,2* ion from 
methane is relative to the CH,* ion; 
CH.+ and CHD? ions from mono- 
deutero-methane are relative to the 
CH;D* ion. Ionizing current is con- 
stant. 
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20 32 24 26 28 30 
IONIZING VOLTAGE 


does exist an activating influence on H, assuming 
the absence of C2H, in the sample. 

The spectrum for mono-deutero-ethane has 
been tentatively examined and similar effects are 
indicated. i 

In addition to the spectra of CH;D and CH, 
for 55-volt electrons, spectra for electron energies 
ranging from 14 to 32 volts have been obtained. 
™p was calculated for each of these voltages, 
with results as shown in Table VIII. Values for 
tp below 18 volts are not precise since the ion 
currents at these low voltages are small and 
error is correspondingly greater. The pairs of 
Tp values at 17, 16, and 15 volts are those 
calculated from two spectra. 

™p remains remarkably constant (0.40+0.01) 
for the wide range of electron impact energies 
between 55 and 22 volts; at lower energies the 
value increases slightly. It is below about 23 
volts that the ions CH*+ and CD* are no longer 
observed, as may be seen in Fig. 3. This suggests 
that our assumption of the constancy of rp for 
formation of all types of ions is not completely 
justified, and that zp is somewhat higher for the 
heavier ions. 

Ion intensities of individual ions from CH;D 
have been calculated from mp for each voltage 
with the result that ionizing efficiency curves 
may be drawn for each ion. In Fig. 1 the curves 
for the ions (CH,D*) and (CH;+) from CH;D 
and (CH;*) from CH, are plotted. Figures 2 and 
3 are similar curves for the (CX.2*) and (CX*) 
type, respectively. 

The double inflection of the curve for (CH2*) 
from CH, in Fig. 2 is caused by the fact that the 
(CH2*) ion is formed in two different ways, as 


32 34 36 38 


shown by Smith.* The appearance potentials of 
the ions from CH, agreed within experimental 
error with the values given by Smith. The curves 
show that within the experimental uncertainty 
of +0.5 volts the corresponding ions of the two 
molecules have the same appearance potentials. 

The existence of these effects suggest that 
caution must be used in analyzing for deuterium 
isotope concentration with the mass spectrom- 
eter. In some cases it may be necessary to 
obtain the spectrum for the pure tagged com- 
pound before the analysis of a mixture of tagged 
and untagged molecules can be performed. For 
example, if the predicted pattern of Table II is 
used instead of the correct pattern of Table | 
for prediction of concentration of CH;D in a 
mixture of CH;D and CHg, considerable error 
may be made. The analysis for CH3D in a 50-50 
mixture will be in error by 14 percent of the 
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Fic. 3. Voltage of the ionizing electron versus relative 
ion intensities of CX*+ type ions. The CH* ion from 
methane is relative to the CH,* ion; CH* and CD* ions 
from mono-deutero-methane are relative to the CH3D* ion. 
Ionizing current is constant. 


$L. G. Smith, Phys. Rev. 51, 263 (1937). 
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CH;D present. However, parent ions (i.e., those 
formed by the removal of an electron from the 
parent molecule) of corresponding molecules 
have shown the same intensity per unit pressure 
in all cases observed. Therefore, if it is possible 
to use only the parent ions as measures of 
concentration, it should not be necessary to 
examine the spectrum of the pure isotope- 
substituted compound. 

No detailed explanation has been developed 
for these effects. However, the 7 and I deduced 
here give definite information about a specific 
type of reacting system. It seems reasonable to 
suppose that they can be of considerable assist- 


ance in extending present methods of making 
rate and reaction mechanism calculations. 


Note added in proof: Since this manuscript was submitted 
for publication, it has come to our attention (through Dr. 
D. P. Stevenson) that Bleakney, Condon, and Smith, 
J. Phys. Chem. 41, 197 (1937), account for the difference in 
ratio of atomic to molecular ions between Hz and D; on the 
basis of the difference in zero point amplitudes of vibration 
and the Franck-Condon principle. This doubtless provides 
the basis for the explanation of the effects reported above. 
It seems certain that if sufficiently accurate information 
were available on ground normal vibrational wave func- 
tions and on excited electronic states it would be possible 
to calculate the w and the T effect for CH;D and other 
deuterium containing molecules. 
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The mercury sensitized reaction between hydrogen and 
oxygen in a flowing system was studied. The quantum 
yield of the reaction was found to be less than unity, and 
no evidence of a chain reaction at 40° was found. Hydrogen 
peroxide was the principal product, and yields of peroxide 
above 90 percent were found for many of the experiments. 
It was shown that in the gas phase reaction between ozone 
and hydrogen peroxide, water was rapidly formed. The 
yield of product was found to increase with increasing 
oxygen flow rate, with increasing total flow rate for a fixed 
hydrogen to oxygen ratio, and to increase and then decrease 
with increasing hydrogen to oxygen ratio. The reaction 
scheme proposed is 


Hg+hv—Hg"*, 


INTRODUCTION 


T has been shown by Coehn and Grote’? that 
hydrogen peroxide could be detected among 
the reaction products if a mixture of hydrogen 
and oxygen was passed rapidly through an irra- 
diated quartz vessel. In a static experiment or 


1 The experimental work reported herein was carried out 
at the University of Illinois, and this report is based on 
work performed for the O.S.R.D. under Contract No. 
OEM-sr-1452. 

2G. K. Rollefson and M. Burton, Photochemistry and 
the Mechanism of Chemical Reactions (Prentice-Hall, Inc., 
New York, 1939), p. 173. 


Hg*+H.—>H+H+Hg, 
Hg*+0.—02+Hg, 
H+0.—H0Oz, 
H+M—HM, 
HM+HM—H:2+M, 
HO.+HO.—H:02 +0Ox>. 


This leads to the rate equation, 
d(H202) __—s/ka’’ (Oz) Taki /k2(H2) 
dt —hes/ha’”(Oz)-+1 bi /ko(Hz) + (O2)’ 
using the flow rates of hydrogen and oxygen. k:/k: is the 
ratio of quenching efficiencies of hydrogen and oxygen for 


activated mercury or 1.62. k3/k,’’ is taken equal to 
1.2 10-3, 








when using low rates of flow no peroxide was 
detected. That hydrogen and oxygen can react 
in the presence of photo-activated mercury 
atoms at room temperature was first established 
by Dickinson.* Later it was demonstrated by 
Marshall‘ that, in a flowing system, hydrogen 
peroxide was formed and that the quantum 
yield of the reaction was about 6.6. This value 
was later corrected to about 2.5.5 Bates and 

3 R. G. Dickinson, Proc. Nat. Acad. Sci. 10, 409 (1924). 


4 A. L. Marshall, J. Phys. Chem. 30, 1078 (1926). 
5 A. L. Marshall, J. Am. Chem. Soc. 54, 4460 (1932). 
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Fic. 1. Gas flow system. 


Taylor® and Bates and Salley’? demonstrated that 
under suitable flow conditions, and wit a proper 
ratio of hydrogen and oxygen, hydrogen peroxide 
was the principal reaction product. Franken- 
burger and Klinkhardt® report results at variance 
with those of Marshall and Bates and co-workers 
in that they found a quantum yield of only about 
unity and considerable water in the reaction 
products. 

The data existing in 1933 were considered 
critically by Bates® who advanced the following 
reaction scheme: 


Hg+hv—Hg* I, (1) 
He*+H.—-H+H+Hg ky (2) 
Hg*+0,—02+ Hg ke (3) 

H+0,—-HO; ks (4) 


HO.+H.—H,0.+H kg (5) 
HO.+HO.—H:02+ Osx. ks (6) 


This chain mechanism allows quantum yields of 
greater than unity and predicts that hydrogen 
peroxide is the sole product. The mechanism is 
based on assumptions in regard to the incident 
intensity in the experiments of the earlier inves- 
tigations, and explains the difference in the 


6 J. R. Bates and H. S. Taylor, J. Am. Chem. Soc. 49, 
2438 (1927). 

7]. R. Bates and D. J. Salley, J. Am. Chem. Soc. 55, 
110 (1933). 

8 W. Frankenburger and H. Klinkhardt, Zeits. f. physik. 
Chemie B15, 421 (1932). 

9]. R. Bates, J. Chem. Phys. 1, 457 (1933). 
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quantum yields reported by Frankenburger and 
Klinkhardt, and by Marshall quantitatively by 
use of the differences in incident intensity. The 
quantum yield given by Frankenburger and 
Klinkhardt is based on the use of monochloro- 
acetic acid as an actinometer and the values of 
the. quantum yield for the photo-hydrolysis 
reaction given by Rudberg.!® However, more 
recent work" has established a value somewhat 
less than that reported by Rudberg. These newer 
data indicate that the values given by Franken- 
burger and Klinkhardt should be lowered, prob- 
ably well below unity. These considerations, 
especially the uncertainties in quantum yield, 
throw some question on the necessity for pos- 
tulating a chain mechanism. In this paper a 
re-investigation of some of the experimental 
work and a revised interpretation of the reaction 
mechanism are presented. 


EXPERIMENTAL 
Apparatus 


The apparatus used to study the photo- 
chemical reaction is shown diagrammatically in 
Fig. 1. Tank oxygen and hydrogen were measured 
by calibrated flowmeters, F; and Fy». The hy- 
drogen was made oxygen free by prior passage 
over heated copper. Constant flow was main- 
tained by allowing a small amount of gas to 
escape through the constant head stabilizers, S; 
and Ss. The gases were mixed and dried in D 
over magnesium perchlorate. The mixture of 
gases then passed through the mercury saturator 
M-—C. From here they were brought to 20° in 
the water bath, C, and then passed through the 
quartz reaction tube R at 40°. The reaction tube 
was enclosed by a mercury resonance lamp of the 
rare gas-mercury discharge type which was in 
the form of a helix 25 mm in diameter and 200 
mm long, L. Condensable product was collected 
in the liquid air trap N. 

To insure a constant pressure of mercury 
vapor it was decided to saturate the gas stream 
at some temperature below the reaction tem- 
perature so that there would be no condensation 
of mercury vapor in the reaction tube. Mercury 
is held in cuplike depressions, 6. The electrical 

10 E. Rudberg, Leitz. f. Physik 24, 247 (1924). 


1 R. N. Smith, P. A. Leighton, and W. G. Leighton, 
J. Am. Chem. Soc. 61, 2299 (1939). 
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heater, c, and the water-jacketed condenser, a, 
maintain boiling water in the space surrounding 
the cups. The gas stream is very likely only 
partially saturated at 100° in passing over the 
mercury containers. However, reducing the tem- 
perature to 20° in C insures saturation at this 
temperature, and a constant pressure of mercury 
vapor is obtained. 

It is difficult to trap out a small quantity of 
condensable gas from a large quantity of non- 
condensable gas at atmospheric pressure. Using 
liquid air or dry ice as the refrigerant, fogs are 
readily formed and are carried through the trap. 
In the case where water is the condensable gas, 
a characteristic snow may sometimes be seen in 
the effluent gas stream. Bates and Salley’ in 
experiments similar to ours found that it was 
necessary to use several U type or bubbler type 
traps in series in order to insure capturing almost 
all of the condensable fraction. In our experi- 
ments it was found that it was not possible to 
collect all of the product in an ordinary bubbler 
type trap, and a trap of rather unique design 
was constructed. A stream of room temperature 
air enters the well as shown. This sets up a large 
temperature gradient in the annular space 
through which the gases are forced to pass. The 
space is only about 1 mm thick. The arrangement 
is effective in preventing or dissipating water or 
hydrogen fogs, and the entire condensable 
product is collected on the cold surface of the 
outer glass cylinder. 


Analytical Procedure 


Although considerable work has been carried 
out to establish the fact that hydrogen peroxide 
could be produced in yields in the neighborhood 
of 85 percent or even higher, the evidence did not 
appear to be conclusive. Bates and Salley have 
described their technique in considerable detail. 
At the end of an experiment, the trap was allowed 
to come to room temperature, weighed, and then 
titrated. Since the gas mixture was usually about 
10 volumes of hydrogen to 1 of oxygen, the dif- 
ference between the weight of this gas mixture 
and the weight of an equal volume of air would 
be an important factor. In this study it was 
found that diffusion was too slow to allow for 
replacement of the gas mixture by air. However, 
dry air could be used to flush out the trap while 


the trap was still very cold so that none of the 
condensed product was lost. The most satis- 
factory technique was to turn off the hydrogen 
gas at the end of an experiment and then con- 
tinue the flow of oxygen until all of the hydrogen 
was replaced. A relatively small correction for 
the difference in density of oxygen and air could 
then be made. This procedure served a dual role 
in that it eliminated the danger of leaving an 
explosive mixture in the reaction system. The 
peroxide content was then determined by 
titration with potassium permanganate. 


Reaction with Mercury 


In many of the experiments, an oxide deposit 
was observable in the region just past the 
illumination zone. In all experiments, however, 
the major part of the mercury was collected in 
the trap along with the condensable products. 
This mercury appeared black until the trap was 
allowed to warm up, when it assumed an ap- 
pearance characteristic of mercuric oxide. Bates 
and Salley have observed similar effects. It was 
important to know whether the reaction resulted 
in catalytic decomposition of the peroxide. 

In experiments in which solutions of 30 percent 
hydrogen peroxide were mixed with mercury, 
mercuric oxide was forned and only stoichio- 
metric amounts of peroxide were used up; 
oxygen gas is liberated in the reaction. Hydrogen 
peroxide was found to decompose only very 
slightly more rapidly over mercuric oxide than 
it does ordinarily in Pyrex vessels. These results 
are similar to those obtained by Elder and 
Rideal’ for the reaction between hydrogen 
peroxide gas and liquid mercury. 

It was still necessary to investigate the pos- 
sibility of catalytic decomposition of hydrogen 
peroxide when high percentage liquid peroxide 


TABLE I. Effect of warming the trap on hydrogen 











peroxide yield. 

Expt. Trap condition Product H202 H202 
No. when rinsed wt., mg wt., mg percent 
74 Room temp. 100 79.7 79.7 
75 Cold 81.6 
76 Room temp. 98 80.3 82.0 
77 Cold 79.9 








2]. W. Elder, Jr. and E. F. Rideal, Trans. Faraday Soc. 
23, 545 (1927). 
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Fic. 2. Yield of product versus oxygen flow rate. 
Hydrogen flow fixed at 3000 ml/min. 


is in contact with liquid mercury. This could 
account for appreciable peroxide loss on warming 
the trap prior to weighing the final product. If 
the reaction results only in a stoichiometric loss 
of hydrogen peroxide caused by the reaction 
with the small quantity of mercury present, the 
amount decomposed should be negligible. Any 
decomposition loss should be avoided by rinsing 
out the trap immediately after removing from 


TABLE II. Variation of O2 flow rate at constant He flow 
of 3000 ml/min. 








flow Prod- H202 H:20 Moles/ 
Expt. ml/ uct H2O2 moles/ moles/ r 
No. min mg % hr hr total #®H:02 7 





115 300 68 94.8 0.0114 0.0012 0.0126 0.27 0.29 
118 300 60 91.2 .0097 .0017 .0114 23 

121 300 64 85.6 .0097 .0030 0127 -23 30 
125 300 73 += 82.0 .0106 .0043 .0149 25 35 


. J J ° 48 
128 2000 104 70.0 .0129 .010 0229 30 54 
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Fic. 3. Yield of product versus total flow rate. 
Hydrogen to oxygen ratio fixed at 3 to 2. 


the liquid air and titrating the hydrogen peroxide. 
Since the trap cannot be weighed, it is not pos- 
sible to determine the percent of peroxide in the 
reaction products when analysis is made in this 
manner. It may be seen from Table I that there 
is no observable loss caused by the interaction 
between mercury and hydrogen peroxide in the 
trap. 


The Quantum Yield 


A determination of the incident quanta was 
made using uranyl oxalate as an actinometer.™ 
The reaction vessel was filled with a solution 


TABLE III. Variation of total flow rate at fixed H2:O: ratio. 
Ratio of H2:02=3:2. 











Total Prod- 
Expt. flow uct H202 H202 H20 Moles/hr 
No. ml/min mg % moles/hr moles/hr total 
139 1000 = 31 56.5 0.00308 0.0043 0.00738 


138 2000 39 = 65.5 .00450 0043 00880 
148 2000 36 8 68.3 00433 0037 -00803 
136 §=6§3300 )3=554—S 70.9 00675 0053 .01205 
149 3300 55 67.8 .00657 .0060 .01257 
143 5000 68 73.8 00893 0060 01493 
146 5000 87 = 80.3 .0123 .0057 .0180 

150 500 83 = 64.7 -00945 0097 .01915 








13 W. G. Leighton and G. S. Forbes, J. Am. Chem. Soc. 
52, 3139 (1930). 
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containing 0.01 mole per liter uranyl oxalate and 
0.05 mole per liter oxalic acid. After a known 
period of irradiation, an aliquot of the solution 
was titrated with standard 0.1 N potassium per- 
manganate. The intensity of light entering the 
reaction cell was found to be 7.110~ Einstein 
per minute. The quantum yields subsequently 
given are based on this value. 


The Hydrogen-Oxygen Reaction 


For experiments in which the flow rate of 
hydrogen was kept constant and the oxygen flow 
rate was varied, results are shown in Table II 
and Fig. 2. 

The yield of product at varying total flow rate 
for a constant ratio of H.O2 was followed. The 
period of irradiation again was 10 minutes. 
Results are shown in Table III and Fig. 3. 

Table IV and Fig. 4 show results obtained for 
a variation of H2:QO, ratio at a constant total 
flow rate of 3300 ml/min. The exposure time for 
these runs was also 10.0 minutes. 


Reaction Between Ozone and 
Hydrogen Peroxide 


Ozone may be formed by the photo-sensitized 
reaction of oxygen with mercury.“ The reaction 
between ozone and hydrogen peroxide has been 
studied in aqueous solution by several inves- 
tigators. The most extensive study was made by 
Taube and Bray,” and the following mechanism 
was advanced : 











H,02+0;—-HO-+ HO:2+0Os2, (7) 

HO2+0;—-HO+20,, (8) 

HO+0;—HO2+0sz, (9) 
TABLE IV. 

Expt. H2O2 Product H2O2 H2O2 H20 Moles/hr 
No. ratio mg 0 moles/hr moles/hr total 
142 10 46 74.0 0.00600 0.0043 0.00103 
137 3 5S faa -00713 .0050 .01213 
136 1.5 54 70.9 .00675 .0053 .01205 
149 1.5 55 6738 -00657 -0060 .01257 
140 0.67 52 # 57.8 -00532 .0073 .01262 
141 0.33 40 25.9 -00184 .0100 .0118 








14 R. G. Dickinson and M. S. Sherill, Proc. Nat. Acad. 
Sci. 12, 175 (1926). 

16H. Taube and W. C. Bray, J. Am. Chem. Soc. 62, 
3357 (1940). 


Reaction (7) is the chain initiating step. (8) 
and (9) result in chain decomposition of ozone. 
(8) and (10) constitute a chain for the reaction 
between ozone and peroxide. 

Since it was likely that ozone was being formed 
in our experiments, we were interested in deter- 
mining whether the reaction occurred in the 
vapor phase and whether it could occur between 
a condensed phase of peroxide and ozone gas at 
low temperatures. Ozone was prepared by 
passing a measured amount of oxygen gas 
through a quartz tube under irradiation of a 
strong mercury discharge. Hydrogen peroxide 
was prepared by the mercury sensitized hydro- 
gen-oxygen reaction. The streams containing 
ozone and peroxide were brought together in a 
mixing chamber, and the resulting product was 
condensed in a dry ice-alcohol trap. The results 
are given in Table VY. Experiment M-1 shows the 
rate of ozone formation in the apparatus, and 
experiments M-3 and M-5 show the rate of 
product yield when no ozone is added. M-4 and 
M-6 show the effect of adding ozone. M-7 shows 
the effect of passing ozone over hydrogen 
peroxide previously trapped out. 

It may be seen that ozone and HO, are present 
at the time of mixing in almost equal molal 
quantities. The reaction between ozone and 
H.O2 results in almost complete removal of 
H.O2 when the two are brought together in the 
gas phase. Passing ozone over H2O» at dry ice 
temperature has no appreciable effect. 
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A very interesting observation is the fact that 
there seems to be a loss of hydrogen in the 
reaction products when hydrogen peroxide and 
ozone are mixed in the gas phase. A comparison 
of M-3, 5, 7 with M-4, 6 shows this to be the 
case. A discrepancy between the total moles of 
peroxide plus water for the experiments with 
and without adding ozone seems greater than can 
be accounted for by experimental error. This 
implies that the reactions involving the free 
radicals formed in the ozone-peroxide reactions 
produce gaseous hydrogen. 


MECHANISM OF THE MERCURY SENSITIZED 
HYDROGEN-OXYGEN REACTION 


The reaction scheme proposed by Bates® leads 
to the kinetic equation, 


d(H2O2) k4(H2) Rila(H2) ’ 
a @® Seen, 
kiJa(H2) 
ki(Hs)+k2(O2) 








This mechanism purports to explain the rate 
variation with the concentration of hydrogen and 
oxygen and the variation of quantum yield with 
intensity of light. Objections to this correlation 
have already been formulated in the introduction 
to this paper. The quantum yields reported 
herein are also well below unity which, although 
not proof of a non-chain mechanism, allows the 
formulation of a non-chain mechanism at the 
temperature used. 

Bonhoeffer and Haber!® and Frankenburger 
and Klinkhardt® arrive at a mechanism which 
substitutes, for reactions (5) and (6) above, the 


TABLE V. The ozone-hydrogen peroxide reaction. 








Moles per hour 104 








Expt : 

No. H:0: H:O 0; Total _— 

M-1 0 0 1.59 1.6 Only ozonizer operating 
M-3 1.64 2.4 — 4.0 Ozonizer dark 

M-4 0.19 2.8 — 3.0 

M-5 1.44 2.0 — 3.4. Ozonizer dark 

M-6 013 30 — 3.1 

M-7 1.39 2.1 — 3.5 Ozonizer operating only 


after H2O2 synthesis 
stopped. 








%K. F. Bonhoeffer and F. Haber, Zeits. f. physik. 
Chemie 137A, 363 (1928). 
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reactions 


HO.+H.—H:,0+0OH, (11) 
OH+OH—H.0,, (12) 


and to explain quantum yields slightly above 
unity, 
HO.+H,0—H,0.+0OH. (13) 


Much of the scheme is derived from the high 
yields of water found in.their experiments. Bates 
and Salley’? have shown that this formulation 
requires at least 51.4 percent water in the 
products. The high percentage yields of peroxide 
found in this and the earlier investigations are a 
refutation of this mechanism. 

Surface effects, the complexity of exact 
kinetics in a flowing system, and the obscureness 
of the role played by activation entering the 
oxygen molecule make it difficult to establish 
the mechanism of the reaction. However, it is 
possible, on the basis of some reasonable assump- 
tions, to arrive at a rate equation which seems 
to have some validity. 

The results of Tables II and III indicate that 
the total moles of product depend upon both 
the flow rate and the composition of the gases. 
The fact that the total number of moles of 
product for a fixed hydrogen to oxygen ratio 
depends on the flow rate shows that hydrogen 
is being reformed by interaction between free 
radicals or H atoms. 

Since it is possible to achieve yields of hy- 
drogen peroxide upwards of 90 percent, we may 
assume that all water found is formed from 
hydrogen peroxide. This water may come from 
peroxide decomposing by an over-all equation, 


2H20.—-2H,0+ Oz, (14) 


or water may come from the reaction between 
hydrogen peroxide and ozone, to which reference 
has already been made. If we assume that these 
water forming reactions are independent ones, 
then we can use the number of moles of product, 
water plus peroxide, as the number of moles of 
hydrogen peroxide formed in the reaction. 

The adsorption of H atoms on quartz or 
Pyrex is strong, and it may be concluded that 
H atoms reaching the wall remain and eventually 
recombine to give hydrogen. We may then sub- 
stitute for Eq. (5) the reactions, 


H+M-—-HM,  k,’ (15) 
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the adsorption of H atoms on the walls, and 
HM+HM—H.2+M, (16) 

the adsorption of H atoms on the walls, and 
HM+HM-—H:2+M, (16) 


the formation of hydrogen on the walls. This 
formulation eliminates the chain character of the 
mechanism of Bates and, using the steady-state 
method, leads to the kinetic equation, 


d(H2O2) k3(Oo) Taki (He) 
dt —_-ks(Ox) +a’M (Hz) +h2(Oz) 





Since the number of collisions of H atoms with 
the walls will be governed by the flow rate, the 
effective wall may then be considered as in- 
versely proportional to the flow rate. For k,’ M, 
the expression k,’’ 1/F, where F is the total flow 
rate, may be substituted. Expressed in terms of 
flow rates, the concentration of each species 
becomes the flow rate at which it is introduced 
into the system divided by the total flow rate. 
In terms of flow rates, the kinetic equation 
becomes, 


d(H2O2) k3/R4’’(Oo) Tak1/R2(He) 
dt‘ Reskea!’(Ox) +1 hei/ko(H2) +(O2) 





k, and k2 are the bimolecular rate constants for 
the quenching of activated mercury by hydrogen 
and oxygen. Zemansky" has shown that the 
relative quenching efficiencies of hydrogen and 
oxygen are 1.1 and 0.68, respectively: kik2in the 
rate expression then is equal to 1.1/0.68 or 1.62. 
I, has been shown to be equal to 0.043 Einstein 
per hour. 

It is hardly to be expected that experimental 
data would be closely fitted by the above ex- 
pression, since variability of walls would be such 
an important factor. However, for data obtained 


17M. W. Zemansky, Phys. Rev. 36, 919 (1930). 
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rate for a fixed hydrogen flow of 3000 ml/min. 


over a small time interval and with consistent 
technique, some success might be expected. 
Using a ratio of k3/k4’/’ =1.2X10-, the rates for 
the condition represented in Table II have been 
calculated. Figure 5 shows a comparison of cal- 
culated rates, the curve, and the experimental 
values, circles. The k3/k,’’ ratio indicates that 
reaction (15) occurs much more readily than 
reaction (3). This is a reasonable result since it 
would be expected that some means of dissipation 
of excess energy is necessary to stabilize HO, - 
formed from H and Os, while adsorption of H 
atoms on the wall should occur on almost every 
collision. . 
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A very thin layer of radioactive copper was electrolytically deposited on a copper blank. 
The surface was then oxidized in air at 1000°C for 18 minutes, giving an oxide layer with a 
thickness of 1.25X10-* cm. After quenching, successive layers of the oxide were removed 
chemically, and the copper activity in each layer was measured. The observed self-diffusion 
of radioactive copper in the oxide agrees quantitatively with a theory based on the following 
assumptions: (a) The oxide grows by diffusion of vacant Cu* sites from the outer surface of 
the oxide inward to the metal. (b) The concentration of vacant sites as the oxygen-oxide 
interface is independent of the oxide thickness, and drops linearly from this constant value to 
zero at the metal boundary. (c) Accompanying the inward flow of vacant sites, there is a flow 
of positive electron holes such as to maintain electrical neutrality. (d) Self-diffusion of copper 
ions takes place only by motion into vacant sites. The results give a fairly direct confirmation 
of the theory of oxidation first suggested by Wagner. 





I. INTRODUCTION 


RENKEL, Schottky, Jost, and particularly 

Wagner,' have explained ionic conductivity 
and diffusion in solids as resulting from certain 
types of lattice defects (departures from perfect 
periodicity) which are in equilibrium with the 
lattice at high temperatures. Wagner? has ex- 
tended these ideas to explain rates of reactions 
in the solid state which are limited by the dif- 
fusion of one of the constituents. 

At high temperatures, copper oxidizes in the 
presence of oxygen gas in the form of a uniform 
film of CusO which, except for very thin films, 
increases in thickness proportionally with the 
square root of the time. The rate of growth 
increases rapidly with increasing temperature, 
and also depends on the oxygen gas pressure.® 
On the basis of a large number of experiments,‘ 


1See, for example, W. Jost, Diffusion und Chemische 
Reaktion in festen Stoffen (T. Steinkopff, Dresden, 1937), 
or N. F. Mott and R. W. Gurney, Electronic Processes in 
Ionic Crystals (Oxford University Press, Oxford, 1940), 
Chaps. II and VIII. 

2C. Wagner, Zeits. f. physik. Chemie B21, 25 (1933); 
B32, 447 (1936); C. Wagner and K. Griinewald, Zeits. f. 
a Chemie B40, 455 (1938). See also C. Wagner, 

rans. Faraday Soc. 34, 851 (1938), N. F. Mott, Trans. 
Faraday Soc. 36, 472 (1940). 

3N. B. Pilling and R. E. Bedworth, J. Inst. Metals 29, 
529 (1923). A logarithmic law is followed at low tempera- 
tures. Cf. B. Lustman and R. F. Mehl, Metals Tech. 8, 
1317 (1941). 

‘The first comprehensive experiments are those of H. 
Diinwald and C. Wagner, Zeits. f. physik. Chemie B22, 
212 (1933). Later experiments, together with a comparison 
of theoretical arid experimental values for the rate of 
oxidation, are summarized by C. Wagner and K. Griine- 
wald, reference 2. 


Wagner and co-workers have developed a theory 
to account for the reaction rate which includes 
the following assumptions: 


(a) The lattice defects in CuO are vacant sites ordinarily 
occupied by Cut ions (so-called Schottky defects). These 
vacant sites have an effective negative charge of one elec- 
tron. 

(b) Electrical neutrality is maintained by an equivalent 
number of positive electron holes (electrons missing from 
the uppermost filled band). These electron holes are re- 
sponsible for the electrical conductivity of Cu,O. 

(c) The electron holes and vacant sites are completely 
dissociated at temperatures of the order of 1000°C. 
There may, however, be incomplete disorder of the Debye- 
Hiickel type. 

(d) The equilibrium number of vacant sites depends on 
the pressure of the oxygen gas in contact with the oxide. 
The concentration in ‘equilibrium with metallic copper is 
small. During oxidation, there is a concentration gradient 
of vacant Cu* sites extending from the oxygen interface to 
the metal. 

(e) The oxide is formed at the oxygen interface. The Cut 
ions arrive by diffusion from the metal through the oxide. 

(f) The mechanism of diffusion of the Cut ions is by 
motion of vacant sites from the oxygen interface to the 
metal. _ 

(g) Electrical neutrality is maintained by an equal 
current of positive electron holes from the oxygen interface 
to the metal. 


The equation for the reaction at the oxygen 
interface may be expressed in the form 


O.—2 Cu 20 +4[Cut],+4[e = Jas 


in which [ ], indicates a vacancy or defect Cut 
ion or electron. Correspondingly, at the metal 
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interface 
Cu+[Cut],+[e—_],—normal lattice. 


The vacant Cut sites are filled by Cut ions from 
the metal and the electron defects are filled by 
electrons from the metal thus eliminating the 
two holes and leaving the normal lattice in their 
place. Except when the oxide layer is very thin, 
the rate of reaction is limited by diffusion of the 
vacant sites from the oxygen interface to the 
metal. A very much simplified version of 
Wagner’s theory is given in Section III. 

It may be appropriate to make some brief 
remarks about the experimental basis for 
Wagner’s assumptions. Wagner and co-workers‘ 
have measured the rate of oxidation, and also 
a number of properties of Cu,O in equilibrium 
with Os, at various pressures, with the following 
conclusions: 


(a) The rate constant for oxidation varies with oxygen 
pressure as p!/7, 
(b) The electrical conductivity® is also proportional to 


pu, 

(c) The stoichiometric excess of O in CuO, as deter- 
mined chemically,* is proportional to p/5, 

(d) The transport number’ for ionic conduction in Cu20 
indicates that vacant Cut sites carry 5X 10~ of the current 
at 1000°C. 


Elementary theory, assuming complete dis- 
sociation of electron defects, and vacant sites, 
predicts that the excess oxygen should vary as 
p’8. Complete association would give a p'/! 
dependence. The observed p'/5 dependence is 
intermediate, suggesting partial dissociation. 
Another possibility, suggested by Wagner and 
Hammer,‘ is that dissociation is complete, but 
that there is incomplete disorder of the Debye- 
Hiickel type. If it is assumed that dissociation is 
complete, so that there are two defect electrons 
for each excess O atom, the number of electrons 
as determined from the stoichiometric excess, 
together with the electrical conductivity, indicate 
a mobility of about 1 cm?/volt sec. at 1000°C.° 
This value is only a little smaller than one would 
predict from extrapolations of mobilities deter- 

5 J. Gundermann, K. Hauffe, and C. Wagner, Zeits. f. 
physik Chem. B37, 148 (1937). 

°C. Wagner and H. Hammer, Zeits. f. physik Chemie 
B40, 197 (1938). 

7J. Gundermann and C. Wagner, Zeits. f. physik 
Chemie B37, 155 (1937). 


* There is apparently a numerical error of a factor of 
ten in the calculations of reference 6. 
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mined from the Hall effect at lower tempera- 
tures. Thus a reasonable fraction of the centers 
must be dissociated. 

As outlined in Section III, Wagner has shown 
that the rate of reaction can be determined from 
the equilibrium number of vacant sites, the elec- 
trical conductivity, and the transport number. 
Assuming complete dissociation, approximate 
agreement with directly measured values was 
obtained. This agreement, together with the 
fact that the rate constant varies with pressure 
in the same way as the electrical conductivity, 
indicates that diffusion during oxidation is 
mainly by dissociated vacant sites. If there is an 
appreciable number of neutral associated centers, 
their mobility must be so small that they do not 
contribute appreciably to the reaction. 

In order to investigate further the diffusion of 
copper ions during oxidation, we have performed 
an experiment using a radioactive copper tracer. 
A thin coating of radioactive copper was electro- 
plated on a copper blank. The surface was 
oxidized in air at atmospheric pressure at a tem- 
perature of 1000°C for 18 minutes, giving an 
oxide film with a thickness of about 0.0127 cm. 
The distribution of radioactive copper in the film 
was then determined. We are here concerned 
with the self-diffusion of copper in CuO which 
has a varying concentration of vacant sites. 

There are two possible mechanisms for self- 
diffusion : 

(a) Motion of copper ions into vacant sites. This 
mechanism will contribute a term to the diffusion coef- 
ficient proportional to the concentration of vacant sites. 

(b) Interchange of two neighboring Cut ions in the 


lattice. The probability should be independent of the 
concentration of vacant sites. 


Our results indicate that the mechanism of 
self-diffusion is (a) rather than (b). Good agree- 
ment between calculated and experimental 
results is obtained if the distribution of vacant 
sites is assumed to vary in position and time in 
accordance with Wagner’s theory of the oxida- 
tion reaction. Our results thus confirm Wagner’s 
theory, without reference to conductivity ex- 
periments. 


II. EXPERIMENTAL PROCEDURE AND RESULTS 


The experimental work described in this paper 
was done in 1938, but publication was delayed 
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TABLE I. Determination of the thicknesses of oxide layers 
from change in weight during etching. 


TABLE II. Percent of radioactive copper in successive 
layers as determined from tests of activity. 











Weight in grams Thickness 
Sample before after diff. in cm 
1 8.097 7.936 -161 is KIO 
2 7.936 7.469 .467 4.67 X 10-3 
3 7.469 7.121 .348 3.48 X 1073 
4* 7.121 6.714 .407 
Oxide in 4** 300 3.00 x 10-3 
Total oxide j a 2.7 Xie 








* This includes some of the underlying copper. 
** By difference. 


because the results are rather incomplete. 
Because of the stress of other work and the 
intervention of the war, further experimental 
work has been postponed indefinitely so that it 
has been decided to publish the available data, 
together with a recent theoretical analysis of the 
problem, at this time.® 

The copper blank used in the experiment was 
electroplated with radioactive copper by J. 
Steigman, then at Columbia University. The 
technique used is exactly the same as that used 
for the preparation of blanks for a study of the 
self-diffusion of copper in the metal, done about 
the same time.!° 

Sufficient radioactive copper (Cu 64, period of 
decay, 12.8 hours) was deposited on one side of 
a copper blank, to give an activity of about 
2000 counts per minute. The blank was then 
oxidized in air at atmospheric pressure at 1000°C 
for 18 minutes, was kept at 500°C for 3 minutes, 
and then quenched in air. The oxide film was 
formed only on the active side of the blank. The 
back of the blank was protected by another sheet 
of copper. At 1000°C the oxide film formed is all 
Cu,0. At 500°C a thin film of black CuO is 
formed on the surface of the Cu,O. The total 
thickness of the oxide was determined by the 
increase in weight of the copper blank after 
oxidation (weight of O2 absorbed). This thickness 
was 12.7 X10-* cm. The CuO layer was less than 
10 percent of this. 

Successive layers of oxide were removed by 
etching the oxidized face of the copper blank in 





® The experimental work was done by the second and 
third named authors. The first named author has been 
concerned only with the theoretical interpretation and 
analysis of results. 

10 The method is described by J. Steigmann, W. Shock- 
ley, and F. C. Nix, Phys. Rev. 56, 13 (1939). See also J. 
Steigmann, Phys. Rev. 53, 771 (1938). 





Counts per Percent 
Wt. of min. per mg of total 
Cu tested Counts corrected for activity 
Sample mg in 5 min. background in sample 
1 21 218 12 22.542 
2 17 180 °1.0 55.0+6 
3 27 155 0.45 17.5+2.5 
4 27 108 0.11 5.0+2.5 








nitric acid. After each etch the blank was washed 
in water, dried, and weighed. These weights and 
the calculated thickness of oxide in each sample 
are given in Table I. The nitric acid and the 
water used in washing the blank after each etch 
were saved. The copper from each successive 
etch was precipitated from the mixture of nitric 
acid and wash water as copper sulfide. Four suc- 
cessive layers were etched off. The last one 
included some of the underlying copper. 

Some of the wet copper sulfide powder from 
each sample was tested for activity by spreading 
it out on filter paper, allowing it to dry, and 
wrapping it around the counter. The results of 
these tests are given in Table II. The total 
activity in each layer is determined by multi- 
plying the activity per mg by the total weight of 
each sample (see Table I). The last column gives 
this activity in percent. The probable errors 
listed were determined from the number of 
counts. Other possible errors are smaller. 


Ill. WAGNER’S THEORY OF OXIDATION 


Wagner’s theory of oxidation, as applied to the 
formation of CuO, may be outlined briefly as 
follows." Let x represent the distance of a point 
in the oxide from the metal surface. Let n,(x) 
and n,(x) represent the concentrations (number 
per cm*) of (dissociated) defect electrons and 
vacant sites, respectively. For an electrostatic 
field F(x) in the oxide, the electronic and ionic 
current densities, from both diffusion and con- 
duction, are 


I,.= —eD,(dn./dx)+o.F(x) (1) 
and 
I,;=eD;,(dn;/dx)+0:F(x), (2) 


11 The following is a much simplified version of Wagner’s 
theory, and is presented to introduce equations needed 
for later reference. Cf. N. F. Mott, reference 2. 
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where D and o represent diffusion coefficients and 
conductivities, respectively. For equal concen- 
trations, D, and o, differ from D; and o; by the 
ratio of the electronic to the ionic mobilities, 
a large factor.” 

Since the space charge is limited, n.~n; 
throughout most of the layer. Each term on the 
right of Eq. (1) therefore differs from the corre- 
sponding term of Eq. (2) by the ratio of the 
mobilities. Furthermore, since there is no net 
current, J,.= —J;, so that J, must be small com- 
pared with either term on the right-hand side of 
Eq. (1). In other words, J, represents a small 
difference between two nearly equal terms. The 
electrostatic field, F(x), of course, adjusts in such 
a way that these conditions are satisfied. Since 
the two terms on the right-hand side of Eq. (2) 
differ from those of Eq. (1) by a constant factor, 
they must also be nearly equal. However, instead 
of subtracting, they add, so that we have ap- 
proximately 


I,;=2eD,(0n;/dx). (3) 


Consequently the effect of the electrostatic field 
in the oxide is to double the current to be 
expected from diffusion alone. Equation (3) is a 
special case of a more general expression derived 
by Wagner, which was in turn based on a 
similar equation for the diffusion of electrolytes. 

Except during the initial stages, the rate of 
oxidation is controlled by diffusion of Cu from 
the metal to the oxygen interface. Diffusion of a 
vacant Cut site from the oxygen interface to the 
metal is equivalent to the motion of one Cut 
ion from the metal to the outside surface. The 
copper ion current (ions per unit area per unit 
time) can be obtained from the current (3) by 
dividing by e, the charge per ion. This current 
must be divergenceless (except for a very small 
contribution used to readjust the value of n; as 
the layer grows). Consequently, d”;/dx is inde- 
pendent of x through the layer. 

Let X(t) be the thickness of the oxide layer 
at time ¢, and suppose that in time dt the thick- 
ness increases by dX. Let N be the normal 
number of Cut ions per unit volume. The number 
of Cut ions which diffuse in time dt is then equal 


12 About 2X 10 at 1000°C (reference 6). 
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to NdX, and making use of Eq. (3), 


dX 2D; On; 


=. (4) 
dt N ox 





It is assumed that the concentrations of vacant 
sites at the metal and oxygen interfaces are 
determined by local reactions and are inde- 
pendent of the thickness of the layer. Denoting 
their difference by An; we have 


on;/dx =An;/X. (5) 


If Eq. (5) be substituted into Eq. (4), we find 
that the rate of growth is inversely proportional 
to the oxide thickness, 


dX /dt=K/X, (6) 
where the rate constant K is 
K=2DAn,/N. (7) 


The thickness at time ¢ as obtained from inte- 
gration of (6) is 


X =(2K#)!=2(DtAn,/N)}. (8) 


The approximate validity of Eqs. (7) and (8) 
has been checked by Wagner and Griinewald.? 
The coefficient of diffusion was obtained from 
Einstein’s relation, eD;=v;kT. The mobility 2; 
was derived from measurements of the transport 
number for ionic conduction, 7;, the electrical 
conductivity, ¢, and the number of vacant sites, 
n; (the latter from the stoichiometric excess of 0) : 


Vi=07;/EN;. 
Thus, in terms of measured quantities, 
K=(2o7r:kT/e?N) (An;/n;). (7a) 


This expression, which is slightly different from 
the one used by Wagner and Griinewald, gives 
good agreement with the observed rate of reac- 
tion. It should be noted that only the relative 
values of the concentration of vacant sites 
appears. An appreciable number of associated 
centers (defect electrons+vacant sites) could 
exist without affecting the agreement provided 
that 


(a) The degree of dissociation does not vary greatly 
with concentration. 

(b) The diffusion constant of associated centers is con- 
siderably less than that of dissociated vacant sites. 








718 

















_— —_— — n 
¥. on 9 
Zz a 
= a | 
a -Tt { 
= 
z | 
us | 
,2) 
4 ! 
8 | 
| 
METAL I Gas 
Xm ' 
x=0 x x=Xi(t) x=X(t,) 


Fic. 1. Schematic diagram of oxide layer showing 
concentration N;(x,?) of vacant Cu* sites at successive 
times. 


IV. THEORY OF DIFFUSION OF 
RADIOACTIVE Cu 


One mechanism by which diffusion of radio- 
active Cu can take place is by motion into a 
vacant lattice site. The vacant site may move 
into any neighboring site, but a Cut ion will 
diffuse only if a neighboring Cut site is vacant. 
This mechanism will give a contribution to the 
self-diffusion coefficient for Cut ions equal to the 
diffusion coefficient for vacant sites, D;, times 
the probability that a given neighboring site is 
vacant, ;/N. This part of the self-diffusion coef- 
ficient may depend both on the position within 
the oxide layer and on the time: 


Din((x, t)/N. (9) 


According to Eq. (5), the concentration of 
vacant sites varies linearly with x, increasing 
from a small value, ”,, at the metal interface to 
Ny =Nm+An;, at the oxygen interface. We may set 


ni(x, t)=nNm+xAn;/X. (10) 


We assume ,, and An; are constant, and that 
the thickness of the layer, X, increases with time 
according to Eq. (8). 

In addition to diffusion by motion into vacant 
sites, it is conceivable that two adjacent Cut 
ions may change places. This would give a 
constant contribution to the self-diffusion coef- 
ficient, indistinguishable from that resulting 
from a uniform distribution of vacant sites, as 
obtained, for example, from the constant term 
mm in Eq. (10). Thus we assume for the self- 
diffusion coefficient, D, for Cut ions: 


; x 


An; 
D=D,+D——. (11) 
N X 
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Replacing X by its value from Eq. (8), 


D=Do+3(Didn;/N)*(x/t"). (12) 


We are now ready to derive the self-diffusion 
equation appropriate during oxidation. Figure 1 
shows the oxide layer at any time ¢. The con- 
centration of vacant sites increases linearly from 
Nm at the metal interface (x=0) to mn, at the 
oxygen interface (x=X). The oxygen interface 
is moving to the right with a velocity 


V=dX /dt=(D,An;/Nt)}. (13) 


The current of radioactive Cut ions at any 
point x in the layer will be the sum of two terms: 
a diffusion current, —D(dc/dx), due to the con- 
centration gradient,-and a drift term, Vc, pro- 
portional to the concentration, which comes 
from the drift of Cut ions from the metal to the 
oxygen interface as the layer grows. Thus the 
total current is 


— D(dc/dx)+ Ve. (14) 


The rate of change of concentration with time is 
equal to the negative of the gradient of the 
current 


0c dD\ dc 
-p—-(v-— —. (15) 
Ox? ox J Ox 
Comparison of (12) and (13) shows that 
V=2(0D/dx), (16) 
so the diffusion equation becomes 
dc d*c OD dc 
—=D—-—-—-_, (17) 
ot Ox? dx Ox 


with D given as a function of x and ¢ by Eq. (11) 
or (12). 

There is a quasi-stationary solution of Eq. (17) 
of the form 


1 ‘ 
cma ete/X), (18) 
which keeps the same shape as the layer increases 
in thickness, with the concentration at similar 
points decreasing inversely with the square root 
of the time. Regardless of the initial conditions 
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existing when the layer is very thin, this solution 
will be approached as the layer increases in 
thickness. This is the solution of interest. 
In order to simplify the resulting equations, 
we set 
s=x/X, (19) 


D,=D,An/N, a=D,/Dj, (21) 


so that we have the following: 





D=D)+D;Z=D,(a+2z), (22) 
Oc 1 1 dg 
—= ——_-g-—_2—, (23) 
at 2t = 2th dz 


0D dc 1Didg 1 dg 


ax ax tt X2%dz 4thdz 
0% 1 1 dg 1 dg 


ax?) X2d2 4D, dz? 





(24) 








(25) 


Substituting into Eq. (17) yields an ordinary dif- 
ferential equation, 

d*g dg 
(z+a)—+ (22—1)—+2g=0. (26) 
dz? dz 
This equation may be expressed in the form 


“( dg 2 1 0 27 
—((e+0)—+2(6- ’e)- Qn) 


The interpretation is that the gradient of the 
current in the relative coordinate system (z coor- 
dinates) is zero. This is the requirement for 
quasi-stationary solution. 

Since the flow of radioactive Cut ions is zero 
at the oxygen interface (g=1), the current itself 
must vanish : 


dg 
(s+a)—+2(z—1)g=0. . (28) 
dz 


In terms of the x coordinate system, Eq. (28) 
implies that 


— D(dc/dx)+ Vce=0 at x=0, 
— D(dc/dx)=0 at x=X. 


There is no flow of radioactive Cut ions from the 
oxide to the metal (x =0), and, of course, no 
flow at the oxide-oxygen interface. 
Equation (28) may be integrated directly to 
give: 
g=A(z+a)* tHe, (29) 
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where A is a constant of integration, determined 
by the total amount of radioactive copper present 
in the oxide. 

In our case, only a negligible fraction of radio- 
active copper will diffuse into the metal; prac- 
tically all will remain in the oxide."* The concen- 
tration would, however, follow the same law 
(Eq. (29)) even if diffusion into the metal were 
appreciable. The constant A would then be 
determined by the fraction of the total amount 
of radioactive copper in the oxide, which would 
depend on the diffusion into the metal, as well as 
on the total amount initially present. Once the 
quasi-stationary solution is attained, the frac- 
tions in the oxide and metal will remain fixed, as 
there is no further flow from the oxide to the 
metal. 

The final solution may be expressed explicitly 
as follows: 


c(x,t) =—{ —+ 


Asx Do\22rtPolrr 
~) e22/X (30) 
UwYXX D, 
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RELATIVE DISTANCE FROM METAL 


Fic. 2. Calculated concentrations of radioactive copper 
in oxide layer based on the following assumptions: Curve 
(A)—(1) self-diffusion takes place only by motion into 
vacant Cu* ion sites, (2) concentration of vacant sites 
increases linearly from 0 at the metal interface to a 
constant value at the oxygen interface (a =0), (3) complete 
dissociation of centers composed of vacant Cu?* sites and 
defect electrons; Curve (B)—same as (A) except that 
centers are undissociated; Curve (C)—same as (A) except 
that there is an additional mechanism of diffusion, giving 
a contribution to the diffusion coefficient independent of 
position in the oxide, and equal in magnitude to that 
caused by motion into vacant sites at the oxide-oxygen 
interface (a=1). 


13 The self-diffusion coefficient of Cu in the metal is of 
the order 10 at 1000°C, as compared with a self-diffusion 
coefficient of the order 10~’, for Cu in Cu2O in equilibrium 
with air at atmospheric pressure and also at 1000°C. 
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TaBLE III. Calculated and observed values of percent 
of total activity in successive layers of oxide. z =fractional 
distance from metal to oxygen interface. D=D,(a+z2) 
=self-diffusion coefficient for Cut ions. Calculated values 
of activity are given for different values of a and for 
assumptions of both complete dissociation and no dissocia- 
tion of centers composed of vacant Cut ion sites and 
defect electrons. 











Calculated 
Range Dissociated Undissoc. 
Sample ofz Observed a=0 0.1 0.25 a= 
1  0.77-1.00 22.5% 21% 20% 18% 17% 18% 
2 0.50-0.77 SS. 54 52% 50% 46% 48 
3 0.22-0.50 17.5 21 22% 23% 25% 26 
4 0-0.22~—s 5. 4 6% 9% 12% 8 








with Dy and D, defined by Eqs. (11) and (20), 
and X, the thickness of the oxide layer given by 
(8). In terms of D,, 


X=2(Dy)!. (31) 


Figure 2 gives a plot of the calculated concen- 
trations of radioactive Cu as a function of the 
relative distance from the metal for three dif- 
ferent assumptions concerning the self-diffusion 
coefficient : 


(A) Complete dissociation of vacant Cu ion sites and 
defect electrons, and Do=a=0. 

(B) No dissociation of vacant Cu ions sites and defect 
electrons, and Dp =a=0. This case is discussed briefly in 
the following section (see Eq. (33)). 

(C) Complete dissociation, but with Do=D,, so that 
a=1, 


It is shown in the following section that (A) is 
in best agreement with experimental data. 


V. COMPARISON OF THEORETICAL AND EXPERI- 
MENTAL VALUES; CONCLUSIONS 


Since the thickness of the samples removed by 
etching were not small compared with the total 
thickness of the oxide, comparison with theory 
is best made by using the fraction of the total 
activity in each sample. If the sample extends 
from z=2, to z=22, the fraction of the total 
activity to be expected can be obtained by 
integration of Eq. (27) 


22 
f (z+ a) 21+@) e-22qz 
21 
f= 





: (32) 
f (sta) 2+a)e—2dz 
0 


For a=0, the integration is elementary. For 
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a0, the integrals can be expressed in terms of 
the incomplete I function, for which tables are 
available. 

Table III gives the experimental values of f, 
expressed in percentages, as obtained from 
Table II, and the theoretical values obtained 
from Eq. (32), for a=0, 0.1, 0.25, and 0.5. 
Complete dissociation of defect electrons and 
vacant sites is assumed. Good agreement (within 
the experimental uncertainties) is obtained for 
a=0. The agreement is almost as good for 
a=0.1, but is definitely worse for larger values 
of a. Thus we conclude that a is most probably 
not greater than about 0.1, and may be con- 
siderably smaller. 

Results of Wagner and Griinewald? indicate 
that the fractional concentration of vacant Cut 
sites in equilibrium with metallic copper is about 
1.2X10-* as compared with about 17Xi0~‘ in 
equilibrium with air at atmospheric pressure, 
both at 1000°C. These values would give: 


1.210 


a=——— = 0.07. 
17X10 


Our results are consistent with this value for a. 

Self-diffusion by interchange of Cut ions in 
neighboring positions must be small compared 
with diffusion via vacant sites, if it exists at all, 
under the conditions of our experiment. 

Some information can also be obtained con- 
cerning the degree of association of defect 
electrons and vacant Cut sites. Without more 
complete data and analysis it is difficult to set 
limits. If there were complete association, the 
centers would migrate by diffusion alone, so the 
factor of two would be missing from Eq. (3). The 
factor of 2 in Eq. (8) would be replaced by v2. 
The factor of 2 would also be missing from Eq. 
(16), and the second term on the right-hand side 
of Eq. (17) would not appear. The analysis can 
carried through for this case, leading to a dis- 
tribution function for radioactive Cu of the form 


g(z) =(z+a)**e~. (33) 

The agreement with observed data is again 
best for a=0, but is not nearly as good as that 
obtained for the assumption of complete dissoci- 
ation. Values obtained from integration of Eq. 
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(33) with a=0 are given in the last column of 
Table III. We conclude that either 


(a) Dissociation is nearly complete or 

(b) Dissociation is not complete, but the mobility of 
dissociated centers is much greater than that of undis- 
sociated centers. 


As far as our results go, they are in complete 
accord with Wagner’s. It is interesting to note 


that the numerical agreement between calculated 
and observed activities is obtained without in- 
troducing any empirical constants and without 
reference to any absolute values for diffusion 
coefficients or conductivities. The shape of the 
concentration curve is determined entirely by the 
mechanisms assumed for the oxidation reaction 
and for the self-diffusion of radioactive Cut. 
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The method used by Gurney and Fowler in explaining the interfacial potential at a metal- 
solution boundary is applied to the zeta-potential. This gives the {-potential as a straight, line 
function of the logarithm of the concentration, in agreement with the shape of the curves as 
found experimentally. To obtain exact agreement between theory and experiment, it is pro- 
posed that the dielectric constant of the electrical double layer is about 0.6 that of the value 
for pure water. This assertion can be checked by accurate measurements on the surface con- 


ductance of solution of electrolytes. 


INTRODUCTION 
SATISFACTORY theory correlating the 


zeta-potential with the concentration of the 
electrolyte in an aqueous solution has not been 
devised. Most of the recent attempts at con- 
structing a theory! have been used on the idea 
of the Gouy double layer. The mathematical 
treatment utilizes solutions of the Poisson- 
Boltzmann equation 


. Are .- ( “<) 1) 
V7 = ——_ )_, 0,2; eX _ ; 
D ‘=1 . kT 





with special boundary conditions. Insufficient 
and unreliable experimental information has 
been largely responsible for the failure of most 
of the attempts. Below, attention is called to a 
method which is believed to be more fruitful 
than previous approaches to the problem. Im- 
provements and modifications of this method, in 

1G. Gouy, J. Phys. [4], 9, 457 (1910); D. L. Chapman, 
Phil. Mag. 25, 475 (1913); O. Stern, Z. Electrochem. 30, 
508 (1924); A. March, Trans. Faraday Soc. 31, 1468 


(1935); A. J. Rutgers and Ed. Verlende, Proc. Kon. Akad. 
v. Wetensch. Amsterdam 42, 71 (1939). 


the light of additional experimental information, 
will probably yield an adequate theory. The 
method employed here has been used by Gurney? 
and Fowler’ in explaining the metallic electrode 
potential, and by Dole‘ in explaining the glass 
electrode potential. 

Recently significant gains have been made in 
regard to the experimental knowledge concern- 
ing the potential existing at the interface 
between vitreous silica and aqueous solutions of 
electrolytes. Hitherto there has been no agree- 
ment among the various experimenters® in this 


2R. W. Gurney, Proc. Roy. Soc. A134, 137 (1931); 
A136, 378 (1932); R. W. Gurney, Ions in Solution (Cam- 
bridge University Press, Cambridge, 1936), pp. 85-88. 

3R. H. Fowler, Proc. Roy. Soc. A136, 391 (1932); 
Trans. Faraday Soc. 28, 368 (1932). 

4M. Dole, J. Chem. Phys. 2, 862 (1934); M. Dole, 
Experimental and Theoretical Electrochemistry (McGraw- 
Hili Book Company, Inc., New York, 1935), pp. 512-515; 
M. Dole, The Glass Electrode (John Wiley and Sons, New 
York, 1941), pp. 261-276. 

5H. Lach and J. Kronman, Bull. Int. de L’Acad. 
Polonaise d. Sci. d. Let. [A], Sci. Math. 289 (1925); H. 
Freundlich and G. Ettish, Zeits. f. physik. Chemie 116, 
401 (1925); H. R. Kruyt and van der Willigen, Kolloid 
Zeits. 45, 307 (1928); Furutani, Kurokochi, and Asoda, 
Jap. J. Gastroent. 2, 148 (1930); H. Lachs and J. Biczyk, 
Zeits. f. physik. Chemie A148, 441 (1930); R. duBois and 
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field, and the reliability of much of the data is 
open to question. Many experimenters have not 
characterized their methods and instruments 
adequately, and, therefore, a valid evaluation of 
their results is difficult. A hopeful indication that 
reliable information has been obtained con- 
cerning these potentials is furnished in the work 
of Rutgers and deSmet® at Ghent (Belgium) and 
Jones, Wood, and Robinson’ at Harvard. This is 
perhaps the first time that systematic agreement 
has been obtained by different investigators, 
using different experimental techniques. Rutgers 
and deSmet, using the method of electro-endos- 
mosis, studied the ¢-potential of solutions of 
metallic nitrates in contact with Jena 16 III 
glass, whereas the Harvard experimenters studied 
solutions of metallic chlorides in contact with 
quartz. The latter group used the streaming 
potential method. It is remarkable that the 
results were practically identical. 

An interesting regularity was observed in both 
sets of data. Instead of the curves passing 
through maxima, as shown in most textbooks,® 
the potentials decreased in a more or less linear 
fashion with the logarithm of the concentration. 
Both groups of workers found out that thorium 
salts did not fit the general behavior pattern. 
Hydrolysis and other effects with thorium salts 
give their solutions properties not common to 
the other electrolytes. In developing a theory, 
one must look for a method of making the zeta- 
potential a straight-line function of the con- 
centration. 


METHOD OF GURNEY AND FOWLER 


The assumptions used in this discussion are 
not entirely those used by Gurney, and the 
mathematical description is largely that em- 


A. H. Roberts, J. Phys. Chem. 40, 543 (1936); A. J. 
Rutgers, Ed. Verlende, and M. Moorkens, Proc. K. Ned. 
Akad. Wetensch. Amsterdam 41, 763 (1938); A. J. Rutgers, 
Trans. Faraday Soc. 36, 69 (1940). 

6A. J. Rutgers and M. deSmet, Trans. Faraday Soc. 
41,-768 (1945). 

7G. Jones and L. A. Wood, J. Chem. Phys. 13, 106 
(1945); L. A. Wood and L. B. Robinson, J. Chem. Phys. 
14, 251 (1946); L. B. Robinson, Ph.D. Thesis, Harvard 
My ao Cambridge 38, Mass., May 1946, p. 125 and 
p. 162. 

8H. A. Abramson, Electrokinetic Phenomena (Chemical 
Catalogue Co., 1934), p. 203; S. Glasstone, Text-Book of 
Physical Chemistry (D. Van Nostrand, 1940), p. 1201; S. 
Glasstone, Introduction to Electrochemistry (D. Van Nos- 
trand, 1942), p. 534. 
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ployed by Fowler. Ions in solution are assumed 
to be affected by electrostatic and thermal forces, 
whereas the ions adsorbed at the surface are 
acted upon by specific adsorption forces in 
addition to the aforementioned types. This con- 
cept has been used before by A. March,’ among 
others. Let Um be the lowest energy level of 
anions adsorbed on the silica surface; it is 
assumed that the phase boundary is impervious 
to cations to the extent that no cations are 
adsorbed on the quartz surface. Let U, be the 
lowest energy level of anions in the bulk of 
solution i.e., energy of the anions without thermal 
motion. The corresponding energy for cations is 
U. but is uninteresting at the moment. The 
anions are able to make transitions from equi- 
librium positions in solution to equilibrium 
positions on the surface of vitreous silica. At the 
same time anions already deposited on the 
surface, having high enough energies, are able 
to make reverse transitions and resume equi- 
librium positions near water molecules. When 
kinetic equilibrium is established, the number of — 
ions depositing on the surface is equal to the 
number leaving the surface dissolving in solu- 
tion (per unit time). It is assumed that the 
Boltzmann distribution prevails both in solution 
and on the quartz surface. 

The problem to be discussed is that of the 
equilibrium between anions in solution and on 
the vitreous silica surface. The anions in solu- 
tion are not entirely free but are bound to water 
molecules or to groups of water molecules. The 
partition function for the anions in solution is 


Octution —_ Ww(T) oP exp [— U;/kT ], (2) 


where W represents the number of water mole- 
cules available for each anion and w(T) repre- 
sents the partition function for each water 
molecule. p; is the statistical weight of each state 
of energy U;; it is quite likely that the per- 
turbation resulting from interacting levels makes 
each p; equal to unity. In the absence of any 
interfacial potential ¢, the partition function for 
an anion on the vitreous silica surface would be 


Qsitica(¢ =0) = obi exp[—U;/kT], (3) 


9 A. March, Trans. Faraday Soc. 31, 1468 (1935). 
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where the symbols have the same meaning as 
before. The existence of an interfacial potential 
¢ shifts each energy level by an amount zef with 
respect to the condition given in (3). Conse- 
quently the correct partition function for the 
ions on the surface is 


Onstea ™ y pyexp[—(Ujtset)/kT]. (4) 


The assumption is also made that the specific 
adsorption forces holding the ions on the surface 
are of the same nature as those which hold the 
units of the vitreous state. Evidence shows that 
the forces existing between particles in the 
vitreous state are similar to those which exist in 
the crystalline state. Hence the adsorbed ions 
may be treated as units of a crystal; the problem 
is now the same as that of a metallic crystal and 
its ions in solution. The additional assumption 
is made that the number of adsorbed anions, on 
the surface of the silica, is a constant independent 
of the concentration of the solution (in the 
limited concentration range in question), i.e. the 
adsorbed ions are fixed to a finite number of 
adsorbing centers. Gurney, Fowler, and Dole 
with their models, assumed that the probability 
of a transition from the surface to the solution 
depended only on the number of water molecules 
in the vicinity of the surface. Here the additional 
assumption is used that the number of transi- 
tions in a given time depends also on the number 
of adsorbed ions. This, however, is a constant 
and contributes nothing to the form of the par- 
tition function. 

Fowler’® shows that the equilibrium condition 
in terms of the partition functions is 


weston 
Qsitica 
Equation (2) may be written as 
Qsotution = Ww(T) (po exp [— Uo/kT ] 
+piexp[—Ui/kT]+---). (6) 
This may be rewritten as 
Qsotution = Ww(T) po exp | — Uo/kT] 


X(1+:1/po exp [—(Ui— Uo)/kRT +--+). (7) 


1 R. H. Fowler, Statistical Mechanics (Cambridge Uni- 
versity Press, Cambridge, 1929), p. 117. 





(5) 


n= 


Another way of writing (7) is 
Qestation = Q'(T) exp [— U./k T], (8) 


where Uy)= U, and Q’(T) contains the rest of the 
expression. Likewise (4) may be written 


Qsities = Q"(T) exp [—(Um+2e5)/kT]. (9) 
Rewriting (5), one obtains 
n Q’(T) exp (— U./kT) 


pre ing ’ (10) 
W Q(T) exp L[—(Un+zef)/kT] 





where /W is, for all practical purposes, the con- 
centration of the anions in solution. 
From (10), it is immediate that 


/ 





) . 
wl Tat + Un— Us) (11) 


Equation (11) is derived by assuming that each 
water molecule is available as a potential 
equilibrium position for an adsorbed anion. 
However, each anion in solution puts out of 
action, say g water molecules, and, consequently, 
the number of water molecules available as 
equilibrium positions is (W—qn), and (11) takes 
a more rigorous form and becomes 


U.—U.. RT c Q(T) 
1—ge Q’ (T) 


2F 2F 

where U, and U,, are molar quantities instead of 
the molecular quantities in (11); R is the molar 
gas constant, and F is the Faraday constant. For 
this type of model, Fowler suggests that 
Q’(T)/Q’'(T) is very nearly unity, but need not 
be for a more general model. Consequently a 
simplification of (12) yields the approximate 
formula for the ¢-potential 


U.—Un RT 
¢=——_—+— Ine. (13) 
2F 


AT} In c—In 





COMPARISON WITH EXPERIMENT 


In calculating the zeta-potential from stream- 
ing potential data, one uses the formula, 


¢=(49n/D)(xE/P), (14) 
and from electro-endosmosis, 
£=(49n/D)(l/E), (15) 
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Fic. 1. The newly defined zeta-potential (1.7 times the 
conventional values) versus concentration. The points 
represent the data of Jones, Wood, and Robinson; the 
solid lines give the slope of the lines in the proposed theory. 
Potassium chloride, ©; barium chloride, @; lanthanum 
chloride, e@. 


where D is the dielectric constant of the electric 
double layer, 7 is the viscosity coefficient of the 
liquid in the double layer, x is the specific con- 
ductance of the solution in the capillary tube, 
l is the length of the capillary, P is the pressure, 
and E is the measured potential. In both cases 
all of the quantities are measured except 7 and 
D. These are, in general, taken as the same as for 
pure water. Some authors!! have pointed out 
that this assumption may be in error, especially 
in connection with the dielectric constant. Con- 
sequently, the values of the zeta-potential 
(Helmholtz-Perrin) as calculated are probably 
erroneous because of the uncertainty in the value 
of the dielectric constant. Thus, exact agreement 
between this theory and the data of Jones, Wood, 
Robinson, Rutgers, and deSmet is not expected. 
The dipole moment per unit area of the electrical 
double layer is the quantity actually measured 
by the experiments, and further assumptions 
must be made in order to convert the moment 


11H. B. Bull and R. A. Gortner, Physics 2, 21 (1932); 
J. W. McBain and M. E. L. McBain, Zeits. f. physik. 
Chemie 161, 279 (1932) ; E. A. Guggenheim, Trans. Faraday 
nce * 5 (1940); L. A. Wood, J. Am. Chem. Soc. 68, 432 

1946). 
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into potential. However, one should expect the 
slopes of the lines representing the experimental 
data to be the same for the various salts. Such 
was actually observed for most of the salts. 

A formula similar to (13) was proposed by 
Rutgers and deSmet in which C referred to the 
concentration of the cation, in contrast to the 
C of (11) and (13) in this article which refers to 
the anion concentration. Their data, with di- 
valent cations, was practically the same as that 
reported by Wood and Robinson.’ The factor 
2.303 RT/zF gives a slope of 29 mv for the 
¢—log C curve in agreement with both groups of 
investigators. If C in (13) referred to the con- 
centration of the cation, then one would expect 
the slope of the curve for the uni-univalent elec- 
trolytes to be greater than any of the others. 
Previously, Rutgers and Verlende! had derived 
the formula, and their last statement in the 
article is precisely: “‘. . . fiir 1-1, 2-1, 3-1, 4-1- 
wertige Elektrolyten sind die ¢—log C Kurven 
also Geraden, deren Neigung bzw. 58, 29, 19.3 
und 14.5 mV. betragt.’”’ But the above assertion 
was not substantiated by either Rutgers and 
co-workers or Jones, Wood, and Robinson. 

This writer believes that the slope of 29 mv 
(actually about 34 mv by Wood and Robinson) 
bears no relation to the factor 2.303 RT/2F. The 
theory proposed in this article makes the slope 
of each line about 58 mv. Therefore, it is pro- 
posed that the data should be multiplied by a 
factor of about 1.7. This proposal means that the 
dielectric constant (or the ratio of D to n) for the 
electric double layer is only 0.6 that of the value 
for pure water. This assertion could be tested 
experimentally by determining the number of 
ions at the surface (not adsorbed on the surface) 
from surface conductance measurements. Veri- 
fication of the theory would result if the number 
of ions calculated in this manner would be about © 
five times the number calculated from the 
Helmholtz-Perrin zeta-potential. 

A graph is given of the newly defined zeta- 
potential (1.7 times the Helmholtz-Perrin zeta- 
potential) versus concentration (see Fig. 1). The 
value for water is greater than that of any of the 
salts.” 


217. A. Wood, J. Am. Chem. Soc. 68, 437 (1946). 
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The absorption of water and dehydration of these films is found to obey a diffusion equation 
by means of which a diffusivity may be defined. This is not a constant but depends on the water 
concentration in the foils. A rapid change in water content immediately after changing the 
ambient conditions suggests the presence of adsorbed surface layers. The total absorption is 
found to increase faster than the relative humidity. 





INTRODUCTION 


OR experiments on scattering of high energy 

protons by hydrogen carried out by the 
authors at Princeton University in 1941 it was 
decided to use pure Cellophane as a scattering 
foil. In order to obtain the absolute value of the 
scattering cross section it was necessary to know 
accurately the percentage of hydrogen by weight 
in the foils. Very pure Cellophane was prepared 
for us by Dr. Nelson Allen of the Cellophane 
Analytical Research Division of DuPont. This 
eliminated the uncertainty due to small amounts 
of glycerine usually added to this material. 
However the Cellophane was found to absorb 
amounts of water depending upon the relative 
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Fic. 1. Weight of }’’ diameter Cellophane foil .0005” 
thick as a function of time it was left in vacuum after 
having been intially in equilibrium with air of relative 
humidity 32.8 percent. 
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*** This work was done in 1940 but not published 
because of war work. 


humidity of the atmosphere and it was necessary 
to understand this phenomenon quantitatively 
in order to know the true hydrogen content of 
the material. 


METHOD AND RESULTS 


Disks of the foil 4’ in diameter (1.27 cm? in 
area) were punched with a special tool from a 
sheet placed between clean paper. These were 
handled with tweezers and weighed in a quartz 
fiber microbalance built along the lines of that 
described by Strong! and put in a chamber 
which could be evacuated. 

The rate of dehydration of the foil in vacuum 
was investigated as well as the rate at which it 
picked up water again when air of a known 
relative humidity was re-admitted to the balance 
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Fic. 2. Weight of Cellophane foil as a function of time 
after having been immersed in air of relative humidity 
34.9 percent. Foil was initially in vacuum. 


1Strong, Procedures in Experimental Physics (Prentice- 
Hall, Inc., New York, 1942), p. 215. 
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Fic. 3. Weight of .0001” Cellophane foils in equilibrium 
with air of various relative humidities. 


chamber. Relative humidities were measured 
with a sling hygrometer. 

Figure 1 shows the weight of a Cellophane foil 
about .0005 inch thick as a function of the time 
it was left in the evacuated chamber after having 
been initially in equilibrium with air at 32.8 
percent relative humidity at approximately 30°C. 
The total loss in weight after 20 hours when the 
weight had become constant was 0.264 mg or 
0.208 mg per square centimeter. 

Figure 2 shows the weight of the Cellophane 
as a function of time after air of relative humidity 
34.9 percent had been admitted to the chamber. 
In this case the total gain in weight was larger 
than the initial loss as expected since the air was 
now of higher relative humidity. 

Figure 3 shows the weight of two different 
samples of Cellophane about 0.0001 inch thick 
as a function of relative humidity of the air with 
which they are in equilibrium. For these thinner 
foils equilibrium was reached in about two hours. 
Only one reading was taken in a day allowing 
sufficient time for the foil to reach a steady state 
at the particular value of the humidity which 
existed. The point at 100 percent was taken by 
placing a dish of water inside the balance 
chamber and sealing it up. In this case the weight 
gain was 0.336 mg. 

Figure 4 shows the weight loss in vacuum of 
polyvinyl alcohol film type B after standing in 
air of relative humidity 57.5 percent and the 
regain of weight when the air was re-admitted. 
Weight change was 0.217 mg. 
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Fic. 4. Weight of }’’ diameter foil of polyvinyl alcohol 
type B during dehydration, and hydration, in air of 
relative humidity 57.5 percent. 


DISCUSSION? 


The diffusion of moisture into and out of the 
foil might be expected to follow the equation 


h?(d?M/dx*) =dM /dt, (1) 


where M is the concentration of water in the foil 
in grams per cm’, x is measured along the thick- 
ness of the foil, ¢ is the time, and A? is the diffu- 
sivity. Since the diameter of the foil is large 
compared to the thickness, it may be treated as 
an infinite slab. The boundary conditions for the 
case of dehydration when the foil whose surfaces 
are in the planes x=0 and x=L, is placed in 
vacuum, are as follows: 








M=M), 0<x<L for t=0, (2) 
M=0, x=OorL, for allt, 
the general solution of Eq. (1),’ 
a i ii (sin yx) 
M=A exp (—/h’7’"t) axe (3) 
then becomes: 
4 ow 1 2r?t(2n+1)?)7 
M=M,- >. exp| | 
T n=0 2n+1 L? p 
(2n+1)) 
Xsin rx = , (4) 





to fit the boundary conditions for dehydration. 


2 For an excellent general discussion of this problem see 
R. M. Barrer, Diffusion in and through Solids (Cambridge 
University Press, New York, 1941). 

3 Ingersoll and Zobel, Mathematical Theory of Heat 
Conduction (Ginn and Company, 1913), p. 105. 
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Fic. 5. Semi-logarithmic plots of weight of Cellophane 
foil during dehydration in vacuum and rehydration in air. 
The diffusivities, h*, are computed from the straight parts 
of the curves at times greater than 50 minutes. The ordin- 
ate of the lower curve is the final weight of water (after 
many hours) minus the weight of water at the time 
indicated by the abscissa. This number is larger than the 
corresponding ordinate for the upper curve because the 
final weight corresponds to 34.9 percent relative humidity 
rather than 32.8 percent which was the initial value in the 
case of dehydration. 


For the case when moist air is again let into 
the chamber around the dry foil, the boundary 
conditions are 


M=M,), x=OorlL, for all ¢, 








M=0, 0<x<L, for t=0. (S) 
The solution then is 
4a 1 
M= Mol 1i--> 
a n=0 2n+1 
h?x?t(2n+1)? (2n+1) 
Xexp |- 7 sin rn} (6) 


For comparison with the experiment we wish 
to know the totai amount of moisture in the foil 
at any time or 


L 
W=area: [ Mdx. (7) 
0 
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Fic. 6. Semi-logarithmic plots of weight of polyvinyl 
alcohol films during dehydration and hydration. The 
diffusivities, h?, are computed from the straight parts of 
the curves at times greater than 100 minutes. The lower 
curve shows the final weight of water (after many hours), 
minus the weight at the time indicated by the abscissa. 


For the case of dehydration this is, integrating 
Eq. (4) term by term, 
8 ow 
W= Wo— > 


a n=0 


Wet alt ; 
{>| le © 


where W, is the initial total weight of moisture. 
For dehydration, integrating Eq. (6) gives 





8 w 


xexp | -{— “(On +1)? i 


For large values of time ¢, terms with »>0 
may be neglected and the diffusivity h? may be 
found from the slope of the experimental plot of 
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log W against time. The terms with ~>0 are 
important only for t<~100 minutes. Thus the 
ratio of the second term to the first term for the 
case of Cellophane for various times is as follows: 


Second term 


T (min.) First term 
0 1/9 
50 0.023 
100 0.005 


The first term is 81 percent of the total sum 
even at t=0. An approximation to the theoretical 
curve may thus be made using only the first 
term of the series for £>100 minutes and in- 
creasing this by 19 percent X(1—#/100) for 
t>~100. This gives too large a value to the 
sum for t< ~100 except at t=0. This correction 
is shown by the dotted lines in Figs. 5 and 6. 

It is seen from the semi-logarithmic plots that 
(1) a single exponential term gives a good fit 
to the data for > ~50 minutes for Cellophane 
and for > ~100 minutes for polyvinyl alcohol; 
(2) increasing the single term by 19 percent 
X(1-—#/100) to correct for all the remaining 
terms at values of t>100 (theoretically too large 
a correction except for =0) does not give values 
on the dehydration curve as large as the experi- 
mental ones at small ¢. 

The second fact noted above has led us to 
conclude (a) the diffusivity h? of water in these 
substances is not a constant but increases as the 
water content increases and (b) there is water 
absorbed on the surface of the foil in layers 100 
or more molecules thick. This can account for 
part of the large deviation of the experimental 
data from straight exponential plots in the case 
of dehydration. Another effect noticed was that 
the weight gain when the foil is in equilibrium 
with air of various relative humidities is not 
linear but increases faster than the relative 
humidity. This is in agreement with the hy- 
pothesis of Peirce‘ of alpha- and beta-water in 


4 Peirce, J. Text. Inst. (T) 20, 133 (1929). 


cotton where alpha-water is chemically combined 
and beta-water is condensed at high relative 
humidity in spaces between the hexose units 
where it is held by capillary forces. 

That the increase .of diffusivity with water 
content cannot be the only effect is shown by 
the hydration curves Figs. 5 and 6 where 
log (Wo—W) has been plotted. At small ¢ these 
curves also are seen to deviate upward from the 
theoretical lines even though in this case the 
water content is increasing with time so the 
diffusivity should initially be small if the above 
conclusion is true. This deviation we believe to 
be caused by surface absorption. To account for 
the observed differences between the extrapo- 
lated dotted lines and the observed weight at 
t=0, about 100 molecular layers would be re- 
quired in the case of polyvinyl alcohol and about 
2000 in the case of Cellophane, assuming the 
true surface area to be just the geometric value 
for a smooth disk. Actually the effective surface 
may be much larger due to irregularities of the 
order of molecular dimensions. Schiff' has shown 
that adsorbed surface layers of the order of 100 
molecules thick can be accounted for by Van der 
Waals’ forces in the special case of helium II. 

The opposite effect, namely initial rapid evapo- 
ration of absorbed surface water in the case of 
dehydration would tend to make those curves 
too steep initially, so the effect of the change of 
h? is less pronounced than at first apparent from 
the initial part of dehydration curves. The real 
change of h? should tend to make the average 
slope of the hydration curves steeper than that 
of the dehydration curves at times from 100 to 
500 minutes, since the foil is dryer in the case of 
dehydration and the average diffusivity will be 
smaller. This effect is apparently smaller than 
the accuracy of the data resulting from the 
subtraction required to obtain the hydration 
curves (Figs. 5 and 6). 


5 L. I. Schiff, Phys. Rev. 59, 839 (1941). 
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In this paper are presented tables of the thermodynamic functions heat content, free energy, 
entropy, and heat capacity of boron trifluoride, boron trichloride, and boron tribromide in the 
ideal gas state at 1 atmosphere from 298.16° to 1000°K. These functions were calculated from 
Raman and infra-red spectra in the case of the two first mentioned and from the Raman 
spectrum of liquid boron tribromide. A satisfactory calorimetric check of the entropy of boron 


trifluoride is presented. 





I, INTRODUCTION 


N recent years much interest has been shown 
in the formation of addition compounds with 
the boron trihalides and the use of the trihalides, 
particularly of boron trifluoride, as catalysts, as 
is witnessed by the reviews of Kastner! and the 
series of papers by G. F. Hennion in the Journal 
of the American Chemical Society. These com- 
pounds have passed from academic to industrial 
interest. 

The spectroscopic and molecular data (Part 
II) are sufficiently complete to permit the calcu- 
lations of the thermodynamic functions of boron 
trifluoride, boron trichloride, and boron tri- 
bromide to the rigid rotator, harmonic oscillator 
approximation with some increase of uncertainty 
in that order (Part III). Low temperature calo- 
rimetric measurements and vapor pressure data 
permit a check of the entropy of boron trifluoride 
(Part IV). 

The following physical constants and chemical 
atomic weights were used in the calculations: 


R=1.98718 cal. deg.' mole, 

hc/k = 1.43847 cm deg., 

N = 6.02283 X 10” mole“, 

Atomic weight of “B= 11.00992, 
Atomic weight of °B = 10.01361, 
Atomic weight of equilibrium B= 10.82. 


The constants are those selected by Wagman, 
Kilpatrick, Taylor, Pitzer, and Rossini.2 The 
chemical atomic weights of boron of masses 10 
and 11 are derived from the physical atomic 
1B. Kastner, Zeits. f. angew. Chemie 54, 273, 296 (1941). 
2D. D. Wagman, ‘J. E. Kilpatrick, W. J. Taylor, K. S. 


Pitzer, and F. D. Rossini, J. Research Nat. Bur. Stand. 
34, 143 (1945). 


weights: “B=11.01292 and “B=10.01631 of 
Bainbridge and Jordan as reported by Livingston 
and Bethe’ and the ratio of the physical to the 
chemical atomic scale, 1.000272. 

Recently Ingraham‘ has satisfactorily deter- 
mined the percentage abundance of the boron 
isotopes. His values of 81.17 percent for "B and 
18.83 percent for B and the chemical atomic 
weights of the isotopes yields 10.822 for equi- 
librium boron as compared to 10.82 of the 
international atomic weights. 


_ Il. SPECTROSCOPIC AND MOLECULAR DATA 


The spectroscopic data for boron trifluoride 
include the infra-red spectra of Bailey, Hale, and 
Thompson® and of Gage and Barker® and the 
Raman spectrum of the vapor by Anderson, 
Lassettre, and Yost7* and by Yost, DeVault, 
Anderson, and Lassettre.7’ Herzberg® has sum- 
marized these data and shown that they are in 
agreement with the plane symmetrical model 
(point group D3). This is also in agreement with 
the zero dipole moment and the electron diffrac- 
tion measurements. The latter by Lévy and 
Brockway’ fix the B—F distance as 1.30+0.02A. 
From the study of the fine structure of the 


3M. S. Livingston and H. A. Bethe, Rev. Mod. Phys. 
9, 245 (1937). 

4M. G. Ingraham, Phys. Rev. 70, 119 (1946). 

5C. R. Bailey, J. B. Hale, and J. W. Thompson, Proc. 
Roy. Soc. (London) A161, 107 (1937). 

6D. M. Gage and E. F. Barker, J. Chem. Phys. 7, 455 
(1939). 

7a T, F. Anderson, E. N. Lassettre, and D. M. Yost, 
J. Chem. Phys. 4, 703 (1936); » D. M. Yost, D. DeVault, 
T. F. Anderson, and E. N. Lassettre, ibid., 6, 424 (1938). 

8G. Herzberg, Infrared and Raman Spectra of Poly- 
atomic Molecules (D. Van Nostrand Company, Inc., New 
York, 1945), pp. 298-9. 

°H. A. Lévy and L. O. Brockway, J. Am. Chem. Soc. 
59, 2085 (1937). 
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TABLE I. Heat content, free energy, entropy, and heat 
capacity of boron trifluoride in the ideal gas state at 1 
atmosphere. Molecular weight = 67.82. o=6. I2=I;=78.8,4 
xX 10-, I, =157.7 10- g cm?. 














H® —H0° —(F° —Ho°) 
‘ing 4 T T So Cp? 
cal deg.~! mole 
298.16 9.338 51.357 60.695 12.063 
300 9.355 51.414 60.770 12.094 
350 9.811 52.892 62.702 12.981 
400 10.257 54.230 64.488 13.764 
450 10.686 55.464 66.150 14.454 
500 11.094 56.611 67.704 15.059 
600 11.841 58.70 70.54 16.046 
700 12.497 60.58 73.07 16.792 
800 13.071 62.29 75.36 17.358 
900 13.572 63.85 77.42 17.792 
1000 14.012 65.31 79.32 18.128 








parallel vg band (Av=0.71 cm), Gage and 
Barker determined this distance to be 1.29A. 
This value leads to J2=J3=78.84xX10-" and 
I,=157.7X10-* g cm’. 

The Raman spectrum of liquid boron tri- 
chloride has been studied by Venkateswaran”; 
Anderson, Lassettre, and Yost’*; and Wagner." 
The latter author confirmed the measurements 
of Anderson, Lassettre, and Yost and with 
greater dispersion found three of the four possible 
lines caused by the chlorine isotopic splitting of: 
the symmetric valence vibration, 71, in approxi- 
mately the separations and intensities to be 
expected. Anderson, Lassettre, and Yost’s as- 
signment of the fundamental frequencies is 
as follows: v1=471; 2v2=902(4%BCls3);  v3(2) 
=946("BCls), =989(°BCl3); and v4(2) =253 
cm,” They also made a more plausible assign- 
ment of the fundamental frequencies for the 
infra-red absorption spectrum of boron trichloride 
vapor, based upon the data of Cassie," as follows : 
v1=471; ve=462("BCl3);  »3(2) =958(4BCls), 
=996(BCls3); v4(2) =243 cm. The dispersion 
of Cassie’s measurements was insufficient to 
permit the determination of the fine structure. 
Lévy and Brockway’s measurement® of .the 
B—Cl distance, 1.73+0.02A, yields I: = I; = 264.3 
X10-“ and J;=528..10-” ¢g cm?. 

10S. Venkateswaran, Ind. J. Phys. 6, 275 (1931). 

41 J. Wagner, Zeits. f. physik. Chemie 193, 55 (1943). 

2 The totally symmetric vibration »;, involving motion 
of the halogen atoms, is independent of the mass of the 
boron isotope. For the doubly degenerate », the isotope 
effect, which should be small, has not been observed in 
the cases of boron trichloride and boron tribromide. 


% A. B. D. Cassie, Proc. Roy. Soc. (London) A148, 87 
(1935). 
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For boron tribromide there are no infra-red 
measurements and the Raman spectrum studies 
have been made with the liquid.’*""* Though 
the later measurements of Wagner give the 
frequencies v;---v4, less than those of Anderson, 
Lassettre, and Yost by 2, 2, 6, and 1 cm™, 
respectively, those of Anderson, Lassettre, and 
Yost have been accepted. They are 11=279, v2 
= 372("BBrs), v3(2) =806(4BBrs), =846(!°BBrs) 
and »4(2) = 151 cm—. Lévy and Brockway’s meas- 
urement of the B—Br distance, 1.87+0.02A, 
yields Iz:=J;=696X10-° and J,;,=1392X10-” 
g cm’. 

According to the Teller-Redlich product rule, 
for such isotopic molecules as "BX; and “BXs, 
the product of the ratios of the frequencies for 
all vibrations of a given symmetry type is equal 
to the square root of the ratio of the reciprocal 
mean masses, i.e., 


v2(BX3) v3v4(°BX3) [A] 
vo(4BX3)  vgvs(BX;) Lu(BX;)1 





The »’s for boron trifluoride are in excellent 
agreement with this rule and v3v4 values agree 
satisfactorily. For boron trichloride and boron 
tribromide for which (u11/u10)? is 1.0441 and 
1.0466, respectively, direct experimental test of 
this rule is not possible, since the v2’s for “BX; 
have not been measured and the »’s have not 
been resolved. The calculated values of the v2’s 
given by Anderson, Lassettre, and Yost’* are 
not quite in accord with these ratios. For the 
present caiculations the Teller-Redlich rule and 
the v; values for “BX; have been accepted and 
values for “BX; have been calculated. The 
former multiplied by the abundance fraction 
0.8117 and the latter by 0.1883 yield for the 
equilibrium BX;: 


VY} v2 v3(2) v4(2) 
BCl;(g) 471 465.8 965.2 243.2 cm 
BBr;(/) 279 375.3 812.0 151.2 


Ill. THERMODYNAMIC FUNCTIONS 


The frequencies given above for boron tri- 
chloride and boron tribromide, and ones calcu- 
lated from Table 75 of reference 8 by combina- 
tion in the relative proportions of 0.8117 and 


14 R, Ananthakrishnan, Proc. Indian Acad. Sci. A4, 74 
(1936). 
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0.1883, have been used with the moments of 
inertia mentioned in Section II to calculate the 
thermodynamic functions of the three trihalides. 
The symmetry number, ¢o, is 6. Nuclear spins 
and the effect of isotopic mixing have been 
neglected. The values of the heat content, free 
energy, entropy, and heat capacity of boron 
trifluoride, boron trichloride, and boron tri- 
bromide in the ideal gas state at 1 atmosphere 
are summarized in Tables I-III. The precision 
of expression of the free energy and entropy 
does not reflect the uncertainty in the B—X 
distances. Thus, e.g., if 73_x were greater by 
0.01A in each case, i.e., equal to the electron 
diffraction value for boron trifluoride, their 
values would be greater by 0.043, 0.034, and 
0.033 cal. deg.—! mole, respectively. 


IV. CALORIMETRIC AND OTHER 
COMPARISONS 


Eucken and Schréder!® measured the heat 
capacities, heats of transition and fusion of boron 
trifluoride from 12° to 173°K, approximately the 
boiling point. The boiling point has been vari- 
ously reported as —99.9 to —101°C.1® Eucken 
and Schréder obtained the entropy of liquid 
boron trifluoride at 154.5°K from their calori- 
metric data. They combined this with a value 
of the entropy of vaporization at 154.5° and the 
saturation pressure, derived from the vapor 
pressure data and equation of Pohland and 


TABLE II. Heat content, free energy, entropy, and heat 
capacity of boron trichloride in the ideal gaseous state at 
one atmosphere. Molecular weight = 117.911. ¢=6. I2=Is 
= 264.3; 10~”, J; =528.6K 10- g cm?. 














H®°—H,° —(F° —Ho°) 
T°K T z S90 Cp® 

cal. deg.~! mole 
298.16 11.262 58.03 69.29 14.97 
300 11.284 58.10 69.38 15.00 
350 11.873 59.88 71.75 15.79 
400 12.415 61.51 73.93 16.43 
450 12.879 62.99 75.87 16.94 
500 13.306 64.37 77.68 17.36 
600 14.04 66.86 80.90 17.99 
700 14.63 69.08 83.71 18.42 
800 15.13 71.06 86.19 18.72 
900 15.54 72.87 88.41 18.94 
1000 15.89 74.52 90.41 19.10 








16 A, Eucken and E. Schréder, Zeits. f. physik. Chemie 
B41, 307 (1938). 

16 Landolt-Bérnstein, Tabellen, third edition, Vol. I, p. 
338; Supplement III, p. 318. 


TABLE IJI. Heat content, free energy, entropy, and 
heat capacity of boron tribromide in the ideal gaseous 
state at 1 atmosphere. Molecular weight = 250.568. o =6. 
I: =I;=696 X10, J; =1392 K 10-” ¢g cm?. 














H® —H° —(F° —Ho®) 
T°K T T So Cp® 

cal. deg.~! mole! 
298.16 12.609 64.88 77.49 16.25 
300 12.631 64.96 77.59 16.28 
350 13.187 66.97 80.16 16.88 
400 13.695 68.75 82.44 17.41 
450 14.130 70.39 84.52 17.80 
500 14.514 71.89 86.41 18.10 
600 15.16 74.60 89.76 18.58 
700 15.67 76.98 92.64 18.88 
800 16.08 79.10 95.18 19.10 
900 16.43 81.01 97.44 19.25 
1000 16.72 82.76 99.47 19.37 








Harlos,!” and with the entropy of compression 
to 1 atmosphere to obtain the entropy in the 
ideal gas state at 1 atmosphere and 154.5°K. 
The entropy change, S°s9s.2—.S°154.5, for the vapor 
was obtained from values of the characteristic 
temperatures of the fundamental vibrations 
supposedly based on the results of Yost, DeVault, 
Anderson, and Lassettre.7” Their result for the 
entropy in the ideal gas state at 1 atmosphere 
and 298.2°K was 61.17 cal. deg.-! mole. 

The calculations of Eucken and Schréder are 
subject to several minor and two major objec- 
tions. Two of their graphically determined areas 
are in error by approximately 0.1 cal. deg. 
mole—!. Inconsistent (and incorrect) values of 
hc/k appear to have been used in obtaining the 
characteristic temperatures. More important is 
the fact that no corrections for non-ideality 
appear to have been made in determining the 
heat of vaporization. The other major objection 
concerns their use of 1.13 cal. deg.—' mole as the 
entropy change involved in going from 142.15° 
to 144.46°K, which purports to be that involved 
in the transition and 


144.46 
f Cp-d ln T. 

142.15 
According to their Table 4 the energy input per 
mole in increasing the temperature 1.128° (in- 
cluding the transition) was 61.0 cal., the contri- 
bution of the normal molal heat was 42.6 cal., 
leaving 18.4 (given as ca. 18) as the true heat of 


17E,. Pohland and W. Harlos, Zeits. f. anorg. allgem. 
Chemie 207, 242 (1932). 
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TABLE IV. Summary of the calculation of the entropy 
of boron trifluoride in the ideal gas state at 298.16°K and 
1 atmosphere, based on calorimetric and vapor pressure 
data to 154.5°K and statistically calculated AS® from 
154.5 to 298.16°K. 








cal. 








deg.~1 

mole! 

Extrapolation, 0-12.6°K 0.55 

Graphical, solid II, 12.6—142.15°K 18.42 

144.45 

Transition and [' Cp(solid I)d In T 0.74 

Fusion, 1014/144.46 7.02 

Graphical, liquid, 144.46-154.5°K 1.65 

Vaporization at 154.5°K, p=0.2185 atmos. 28.55 

Correction for gas imperfection 0.04 
Entropy of ideal gas at 154.5°K and p=0.2185 

atmos. 56.97 

Compression to 1 atmos. — 3.02 





Entropy of ideal gas at 154.5°K and p=1 atmos. 53.95 
5298.16 — S°154.5 6.85 





Entropy of ideal gas at 298.16°K and p=1 atmos. 60.80 








transition. This indicates an average heat ca- 
pacity in this narrow region of 37.8 cal. deg. 
mole. This is in qualitative agreement with 
their statement : “‘In this region the heat capacity 
was abnormally high (the minimum value was 
~35 cal. deg.' mole“).”’ It would appear that 
the entropy change in going from 142.15° to 
144.46°K should be 0.74 cal. deg.—! mole“. 

Kelley'® recalculated the entropy of boron 
trifluoride in essentially the same manner, but 
with satisfactory graphical areas. His result was 
S°o98.1 = 61.4 cal. deg.—! mole. 

Pohland and Harlos'’? determined the vapor 
pressure of liquid boron trifluoride at thirteen 
temperatures in the range 145.9° to 170.3°K. 
Their equation 


log p (in mm Hg) = —1174.4/T+1.75 log T 
—0.0133507 +-8.0536 


is a fair representation of their data, >> (Pobs. 
—Peatca.) for six positive and seven negative 
deviations totaling 5 and 7 mm, respectively. 
Four values in the range 187.5 to 653 mm by 
Fischer and Weidemann" are in good agreement 
with the equation of Pohland and Harlos. Other 
measurements extant in this pressure range’ do 


18K. K. Kelley, U. S. Bur. Mines, Bull. No. 434 (1941). 

19 W. Fischer and W. Weidemann, Zeits. f. anorg. allgem. 
Chemie 213, 106 (1933). 

” QO. Ruff et al., Zeits. f. anorg. allgem. Chemie 206, 59 
(1932); L. LeBoucher and W. Fischer with W. Biltz, ibid., 
207, 61 (1932). 
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not merit equal consideration since the samples 
were not pure or the measurements were made 
with insufficient accuracy. 

If one uses the equation of Pohland and 
Harlos, with the assumption that the vapor 
obeys the approximate Berthelot equation of 
state in the region concerned 
din p 9 pT, pT.? 

ee 
dT 128 p.T | Si 











AH yap = RT? 





=| 
RTI 


where p,=49.2 atmos., ¢,.=—12.25°C,”4 and 
p(liquid) = 2.25—0.003928T, based on Fischer 
and Weidemann’s data, the heat of vaporization 
at the temperature and pressure used by Eucken 
and Schréder and by Kelley is 4411 cal. mole, 
and AS(vap.) =29.55 cal. deg.' mole. This 
replaces 28.68 and 28.74 given by Eucken and 
Schréder and Kelley, respectively. 

The recalculation of the entropy of boron 
trifluoride in the ideal gas state at 298.16°K and 
1 atmosphere is summarized in Table IV. For 
the present calculations the entropy change 
$298.16 —S°154.5 was obtained from the values of 
S°o9s.16 and S°154.5 derived from the spectroscopic 
data, rather than by an integration of calculated 
heat capacities in the region. 

Similar calculations at each of the nine inter- 
mediate temperatures at which Pohland and 
Harlos measured the vapor pressure (i.e., ex- 
cluding the three lowest and one highest temper- 
atures) gave S°29s.15 = 60.69 (average) cal. deg.—! 
mole—!, with six positive deviations totaling 0.37 
and three negative deviations totaling 0.39. The 
agreement with the statistically calculated value 
of Table I, 60.70, is excellent. Eucken and 
Schréder and Kelley’s statistically calculated 


_values were 61.11 and 61.4 cal. deg.-' mole, 


respectively. 

The entropies of boron trichloride and boron 
tribromide calculated from the spectroscopic 
data by Kelley,'® 68.6 and 76.7 cal. deg.—! mole“, 
respectively, are in error as a result of the use of 
the harmonic 272 in place of the fundamental v2. 


21H. S. Booth and J. M. Carter, J. Phys. Chem. 36, 
1359 (1932). 
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The Theory of the Racemization of Optically Active Derivatives of Diphenyl 


eS 
le F. H. WESTHEIMER AND JOSEPH E. MAYER 
Department of Chemistry, University of Chicago, Chicago, Illinois 
id (Received August 2, 1946) 
or 
of T. L. Hill has recently shown that it is possible to compute the magnitude of steric strain in 


certain organic molecules. The present paper summarizes the results of a similar but independent 
investigation directed towards a computation of the rate of racemization of optically active 
(i.e., sterically hindered) derivatives of diphenyl. The energy of the planar form of a sterically 
hindered diphenyl can be approximated by the equation: 


E=2; }aig2+Ai exp (—d:i/p1)+A2 exp (—d2/p2), 
: where the g; are the normal coordinates of the unstrained molecule in question and the exponen- 
| tial terms are approximations (over the limited range of interest) to the steric repulsions of 


the non-bonded groups which repel each other. There are two such exponential terms, for in 
most sterically hindered diphenyls the repulsion is caused by two pairs of ortho substituents. 








od The equation for Eo, the activation energy for racemization, has been found for the symmetrical 
er case by minimizing E; the constants a;, A, and p, and the dimensions of the diphenyl] derivative 
on appear in the expression for Ep. It is also possible to derive equations for the vibration fre- 
quencies of the planar form of the sterically hindered molecule (and, therefore, find the entropy 
ca of activation); these equations are complicated by the essential degeneracy of the vibrations in 
—, question. 
nis 
nd 
I repulsion of the groups which have been pressed 
on . , together. uantitatively, Hill’s equation is 
d N a recent letter to the editor, T. L. Hill' has —— Quantitatively a in Rt 
n ‘ similar to Eq. (1), where E is the energy of the 
. presented the outline of a general method for : : 
‘or -ayyeee ; apa strained molecule in excess of that of a hypo- 
the quantitative (or at least semi-quantitative) : ; : . 4 
ge . - % : thetical molecule, identical in every way with 
evaluation of steric strain.2 The method applies , , : 
of - - the first, except that it has no steric strain. 
i to those molecules where neighboring but non- shi F 
“ d bonded groups have been forced together so that E=)i 3aq?+E(steric). (1) 
they are separated by a distance less than the Here E (steric) is a function of the distance d, 
— sum of their van der Waals radii. (An example between the groups or atoms which are involved 
a d of the type discussed by Hill is the addition in the steric interference. E(steric) may be 
fi compound of boron trifluoride and a, a’-lutidine, negative (attraction) or positive (repulsion) but 
: where the methyl groups in the base are forced for most cases of interest, and certainly in the 
a close to the fluorine atoms in the boron deriva- case studied here, it is positive. The variables q; 
37 tive.2’ Investigations by the same method as are coordinates which for the molecule in ques- 
rh that suggested by Hill have been under way here tion are proportional to the displacements of all 
’ id for the past nine months, with special reference the atoms from their equilibrium positions. 
ore to the rate of racemization of optically active Equation (1) is written without cross-product 
ae diphenyls. It seems therefore appropriate to terms, that is, without terms a;;q:g;. The assump- 
ah supplement Hill’s publication with those results tion is thus made that the coordinates, qi, are 
= obtained here which are of general applicability. mutually independent. This assumption is satis- 
Hill suggested that the energy of a strained fied for the “normal” coordinates of the molecule 
_ molecule is distributed among many bonds and_ but not in general for the increments (decre- 
is by no means concentrated in the energy of ments) to the Cartesian coordinates of the 
le", eR various atoms. The sum in Eq. (1) must extend 
e of 1T. L. Hill, J. Chem. Phys. 14, 465 (1946). See also I. . ; : . 
Dostrovsky, E. D. Hughes, and C. K. Ingold, J. Chem. ~— all those coordinates, gi, changes in which 
1 v2. Soc. 173 (1946). will affect the value of d. The force constant 
2a ° . . ° . 
i ~ 7 Tt. 8 iL ae, ~~ ' 2 edo: associated with each coordinate q; is represented 


J. Am. Chem. Soc. 64, 325 (1942). 


by a;. The values of the coordinates, qi, are to be 
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Br 


Fic. 1. Planar form (the activated complex) of the diphenyl. 


chosen in such a way that E is a minimum. It 
should be noted that the larger the number of 
terms in Eq. (1) the lower will be the minimum 
value of E. This is the usual situation which 
obtains in a variation treatment, and is a neces- 
sary consequence of the fact that the energy 
depends upon a power of the displacement coor- 
dinates higher than the first. 


II 


The rest of this article is concerned with the 
quantitative treatment for the rate of racemi- 
zation of optically active diphenyls such as 
2,2’-dibromo-4, 4’-dicarboxydiphenyl.* The race- 
mization of a compound of this type will occur 
only if the molecule passes through the planar 
configuration (Fig. 1). This will be the ‘‘activated 
complex”’ for the racemization of sterically 
hindered diphenyls. 

It is worth noting in passing that the proposed 
treatment of steric strain assumes that the race- 
mization involves (among other displacements) 
an increase in the angle, ¢, between each ortho 
substituent and the benzene ring. This increase in 
angle will be sterically hindered by a substituent 
in the adjacent meta position. This conclusion is 
consistent with the experimental fact‘ that the 
rate of racemization of some diphenyl derivatives 
is uniformly and considerably decreased by sub- 
stituents in this meta position. Substituents in 
other positions, although they affect the rate of 
racemization, do so to a much smaller extent, 
and they occasionally increase, rather than 
decrease, the racemization rate. 

In the particular example chosen, only the 
hydrogen and bromine nuclei are close enough 


3W. I. Patterson and R. Adams, J. Am. Chem. Soc. 
57, 762 (1935). 

*R. Adams and H. R. Snyder, J. Am. Chem. Soc. 60, 
1411 (1938). 
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together to interfere strongly. (It is assumed that 
the bromine-bromine repulsions are so strong as 
to preclude the possibility of racemization 
through an activated complex in which these 
atoms are adjacent to one another.)The distance 
between a pair of hydrogen and bromine nuclei 
in the activated complex is represented by d; 
the distance between these nuclei which would 
obtain if there were no deviations from normal 
bond lengths and angles is called do. In the case 
under consideration, the energy of the two neigh- 
boring but non-bonded groups is predominately 
repulsive in character. It is convenient to ap- 
proximate this repulsive energy by A exp (—d/p); 
this approximation should be valid for the values 
of d in question. The constants A and p can be 
evaluated by comparing the expression above 
with the energy curve computed from the second 
virial coefficient® for molecules which contain the 
same groups (here hydrogen and bromine) which 
repel each other in the activated complex. 

In the case under discussion (2,2’-dibromo-4, 
4’-dicarboxydiphenyl) each of the two bromine 
atoms must be forced close to a hydrogen atom 
in the activated complex. The repulsive energy 
for any distance, d, will then be twice that cal- 
culated from the energy curve to which reference 
has been made above. The same considerations 
will apply to other cases in which the diphenyl 
derivative is symmetrically substituted in the 
ortho positions. In the case of unsymmetrically 
substituted diphenyls, the sum of the repulsive 
energy for the two different sets of repelling 
groups would have to be used for E(steric) in 
Eq. (1). The energy of activation for the present 
case is that shown in Eq. (2): 


E=)i209?+2A exp (—d/p). (2) 


Here the constants a; are the force constants® 
associated with the normal coordinates qj. 
Normal coordinates are chosen so that both the 
potential and kinetic energies involve only 
square and no cross-product terms. We wish to 
use such coordinates that the potential energy 
involves only squares but we are not primarily 
interested in the kinetic energy terms. It is 

5 See, for example, O. K. Rice, J. Am. Chem. Soc. 63, 3 
“ a discussion of these constants for benzene see E. 
Bernard, C. Manneback, and A. Verlaysen, Ann. Soc. Sci. 


Bruxelles 60, 45 (1940); F. A. Miller and B. L. Crawford, 
J. Chem. Phys. 14, 282 (1946). 
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therefore possible (except in Part III of this 
paper, in which the entropy of activation is dis- 
cussed) to use bastardized coordinates which are 
not true normal coordinates but involve true 
normal coordinates of part of the molecule (the 
benzene rings) and force constants for other 
parts. In fact, the normal coordinates for the 
substituted diphenyls are not now available. 
However, the normal coordinates for benzene 
are known. In order to carry through the cal- 
culations, it will be necessary to assume that the 
force constants for the various deformations in 
each benzene ring are unaffected by the other 
aromatic nucleus. It will further be necessary to 
make some assumption about the magnitude of 
the force constant for the stretching vibration of 
the central carbon to carbon bond in the diphenyl 
derivative. The value of the activation energy 
computed on the basis of these assumptions 
should be the same as that which would be 
obtained from the normal coordinates of diphenyl 
provided only that a correct value for the 
stretching force constant of the central bond has 
been chosen. It should be noted, however, that 
the summation in (2) extends over all the normal 
coordinates in both benzene rings. 

Preliminary calculations have shown that for 
diphenyls which racemize at measurable rates, 
only small changes in bond lengths and angles 
are necessary to attain the planar configuration. 
Because the displacement of each atom from its 
equilibrium position is small, it is possible as a 
good approximation to represent d as a linear 
function of the Cartesian displacement coor- 
dinates, x;. Here, however, it will be necessary to 





set up separate functions for the distances d, and 
d, between the two hydrogen-bromine pairs, 
since the distances between the atoms in the two 
pairs are affected differently by changes in the 
various Cartesian displacement coordinates. 


dy=dotDilixi, de=dotQVil,’xi. (3) 


There is a linear relationship’ (Eq. (4)) between 
the Cartesian displacement coordinates of each 
atom and the normal coordinates, qi. 


j= Loi CisQi- (4) 


(The normal coordinates of the benzene mole- 
cule, for example, are known in terms of the 
Cartesian displacement coordinates of each 
atom.*) Because of the considerations cited 
above, Eq. (3) may be written in the form (5) 


d,=dot+)>d: digi, dz=do+)> i bi'qi, (S) 
where 


b:=Dilicx and b/=>D hk’ cu. (6) 


The sum in Eq. (2) extends over all the normal 
coordinates, g;. There will be some of these coor- 
dinates, however, changes in which do not change 
the internuclear distance, d; for example, this 
applies as a first approximation to the non-planar 
vibrations of the aromatic rings in a diphenyl 
derivative. For such a coordinate, g,, the cor- 
responding constant }, is zero. It will be seen 
below (Eq. (10)) that when any particular },, 
say 0;, is zero, the value of gq, in the activated 
complex is also zero ; the term in Eq. (2) involving 
this coordinate therefore contributes nothing to 
the energy, E, and may be omitted from the cal- 
culation if this for any reason seems desirable. 


The combination of Eq. (2) and Eq. (5) leads to (7) 
E= Di 30g? + Ed exp [—(1/p) 2: big: J+ Ea exp [—(1/p) 2: 5,9], (7) 


where Eay=A exp [(—do/p) ]. Here Edo is the repulsive potential which would have obtained if no 
bond lengths or angles had been distorted. It is now possible to find the minimum value for E; the 
values of the normal coordinates, q;, in the activated state are represented by gq;° 


(0E/dq;) =a3q;° — (b;/p) Edo exp [—(1/p) i bigi® ]—(0;'/p)Eao exp [—(1/p) 2: bgi°J]=0. (8) 


Let 


1 1 
Ea exp | -- a ba0| = Ea exp | -- oe 6140 = Z. (9) 
p p 


Since only definite values of the normal coordinates, corresponding to the ‘‘equilibrium”’ state of the 


7Cf. G. Herzberg, Infrared and Raman Spectra (D. Van Nostrand Company, New York, 1945). 


8 E, B. Wilson, Phys. Rev. 45, 706 (1934). 
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activated complex, occur in Z, it follows that Z itself is a constant for any particular optically 
active diphenyl. Further, from what has already been said about the constants 6 and b’ and about 
the symmetry of the problem, it is clear that the two sums in Eq. (9) must be equal; this will also 
later be apparent from physical considerations. Then 


(6;+b/)Z 2B;Z 
qi= _ ’ (10) 
pa; pa; 





where Eq. (10) defines B;. Z can be obtained in terms of known constants by the numerical solution 
of Eq. (11), which equation results from the combination of (9) and (10). 


Z=Ed exp —[(2Z/p?) >: B?/a;]. (11) 


Substitution of the values of g;° (Eq. (10)) into Eq. (7) yields Eo, the minimum value of E; Ep is, 
therefore, the activation energy for racemization. 
22? B? : 
Eo=>: 20,(gi')?+2Z =—— 2, —— +22. (12) 
pe * as 


Z is therefore analogous to E(steric) of Eq. 1. It is the actual energy of repulsion of one pair of 
neighboring non-bonded groups in the activated state, and is of course much less than Edo. It is 
also clear that, for symmetrically substituted diphenyls, Z must be the same for the two halves of 
the molecule (i.e., Eq. (9) is valid). 

Equation (12) for the energy of activation for the racemization of diphenyl derivatives is not 
precisely correct, since the resonance energy of the molecule has been neglected. However, the 
resonance energy for dipheny] itself, in excess of that for two benzene rings, amounts only to a few 
kilocalories®; the precise value of this energy term is, however, uncertain. When a satisfactory 
estimate of this resonance energy is available, the true activation energy can be found by subtracting 
the resonance energy from Ep. In the calculation of this latter quantity, the observed length of the 
carbon-carbon bond between the benzene rings in diphenyl should be used; this length is needed to 
evaluate the constants /; and /,’, of Eq. (3). 

However, until an accurate estimate of the resonance energy is available, an approximate method 
of calculation must be employed. The resonance in the planar form of diphenyl] decreases the lengths 
of the carbon-carbon bond between the benzene rings. It follows, therefore, that the repulsion 
between the ortho substituents in the planar form of a sterically hindered dipheny] will be increased, 
because of resonance, over that which would have obtained in the absence of resonance.!° This addi- 
tional steric hindrance will diminish the absolute magnitude of the contribution of resonance to the 
energy of the planar configuration. Since the contribution of resonance even before it is diminished 
by steric effects, is small, the net effect on Ey can probably be neglected." It should be noted, however, 
that in any computation where resonance is neglected, the normal (i.e., aliphatic) carbon-carbon 
bond length should be used for the central bond in a diphenyl derivative. 

Approximate numerical calculation of the activation energy for the racemization of 2,2’-dibromo- 
4,4’-dicarboxydipheny] has led to a value of the order of 20 kcal.; a more refined calculation is in 
progress. 


®G. W. Wheland, The Theory of Resonance (John Wiley and Sons, Inc., New York, 1944). 

1 M. Calvin, J. Org. Chem. 4, 256 (1939). 

1! The net effect, however, of resonance is to lower the energy of activation, Eo, although the amount by which Eo 
is lowered is less than the resonance energy of diphenyl (in excess of that in the benzene rings). This conclusion follows 
because (a) resonance will be more important in the planar than in the non-planar structure of any diphenyl and (b) the 
energy of a resonance hybrid is always less (although perhaps only slightly less) than that predicted for the most stable 
individual structure which contributes to the ground state of the molecule in question. The stabilization of the activated 
complex by resonance, therefore, appears inevitable. The opposite viewpoint, which has previously (reference 10) been 
expressed, would appear to be in error. 
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III 
The free energy of activation for racemization is given by the equation 
AF* = Ey— RT In (f*/f), (13) 


where f* is the partition function for the planar activated state and f the partition function for the 
non-planar normal state. The partition function for the normal state can be approximated by the 
product of partition functions, one for each vibrational and one for each rotational degree of 
freedom in the molecule; a similar statement applies to the activated complex. The partition func- 
tions depend, of course, on the vibration frequencies »,* (of the activated complex) and »; (of the 
normal molecule). As a rough first approximation, each term in f* will cancel the corresponding one 
in f. This approximate cancellation is the result in part of the fact that the vibration frequencies in 
the activated complex will not differ very much from those in the unstrained molecule. The can- 
cellation is also the result in part of the fact that, for most of the vibrations in question, hy; is large 
compared with kT; the contribution of such a vibration to the entropy of activation will be small. 
An exception to the approximate cancellation of corresponding terms in the twisting vibration of 
the unstrained molecule about the central carbon-to-carbon bond in diphenyl, since this is the motion 
which, in molecules of high energy, leads to racemization. 

Corresponding terms, however, do not precisely cancel since steric hindrance somewhat modifies 
the vibrational frequencies and energies on the activated complex. The energy of the activated 
complex is in fact represented by Eq. (14), 


E=))i 3a:(9¢°+9,')?+Ea exp [—(1/p) Di bi(gi2+¢.’) ]+ Ea exp [—(1/p) Xi b/ (9+) ] 


(14) 
=>: 24:(¢°+qi’)?+Z exp [- (1/p) di biqi’ ]+Z exp [- (1/p) di bi'q.'], 


where the g;’ are small displacements from the ‘‘equilibrium”’ configuration of the activated complex, 
denoted by q;°. For sufficiently small values of the displacements, q;’, the exponential terms in (14) 
can be expanded in terms of the new coordinates, g,’; this expansion together with Eqs. (10) and 
(12) lead to (15). 


E—Eo= 2s 20:(9i')?+(Z/2p?) Di Db; +-b,'b/)q,'a/, (15) 





where the double summation is carried out over just as the vibration frequencies of the un- 
all values of 7 and’ j, including 7=7. Whenever strained (non-planar) molecule are related to 
the potential energy, E—E», can be represented the corresponding force constants a; by the 
by a quadratic function of this sort, it is in expression 

principle’ possible to eliminate the cross product a;=4n'p?. (17’) 
terms q,’qg;. Further, if the coordinates g; are 
true, normal coordinates, and if they have been 
normalized (i.e., if T= >>; g7), then the elimina- 
tion of cross-product terms from Eq. (15) will 
lead to a new set of normal coordinates such that 


The general transformation of (15) to (16) 
demands the solution of a secular equation with 
n roots, where » is the number of vibrational 
degrees of freedom of the molecule in question. 
The solution of such an equation, in the general 


E-E,y=dia0?2 (16) case, is very difficult, and in many cases cannot 
and reasonably be performed. However, in the event 
T=); Q?. (16’) that the cross-product terms in (15) are small 


compared to the terms a,(q,’)’, the secular 
equation can be solved by perturbation methods. 
The actual solution leads to expressions for the 
force constants a;, the frequencies of vibration 
a= 4n? py (17) »,*, of the activated state, and finally by way 


The frequencies of vibration for the activated 
complex are related to the new force constants 
a;, by the expression 
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of Eq. (13) to the free energy of activation and 
the absolute rate of reaction. Unfortunately, the 
perturbation method is considerably complicated 
by the essential degeneracy of many of the 
benzenoid vibrations. Since for the cases under 
consideration the two benzene rings are iden- 
tically substituted, every vibration of the 
diphenyl molecule is at least doubly degenerate 
and those corresponding to the doubly degenerate 
levels of the benzene ring are fourfold degenerate. 
In principle, these degeneracies can be eliminated, 
and the perturbation method can be carried 
through to completion. However, since the steps 
necessary to this solution are concerned with the 
details of the particular example chosen, it seems 
best to postpone the presentation of the explicit 
solution of this case. The authors hope soon to 


CALDIROLA 


carry out numerical calculations based on this 
development. 


IV 


Although the present treatment has been re- 
stricted to a consideration of the racemization 
of diphenyls, it is clearly of much broader appli- 
cation. An equation analogous to (12) will predict 
the energy of reaction for those cases in which an 
equilibrium constant rather than a rate constant 
can be observed. 


V 
The authors wish to thank Dr. G. W. Wheland 


for his valuable advice and assistance in this 
research. 
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On the basis of the thermo-hydrodynamic theory of the explosive detonation—developed by 
Chapman, Jouguet, and R. Becker—a procedure is developed permitting the deduction of 
the equation of state for gaseous mixtures from explosive processes. The equation is estab- 
lished by starting from several sets of experimental data. In this way one can semi-empirically 
derive the equation of state for gases compressed up to 10*-10° atmos. The results of numerical 


calculations are reported. 


1 


IRECT experimental measurements allow 
one to derive the equation of state for gases 
up to pressures of some thousand atmospheres. 
But we have little knowledge about the properties 
of gases compressed up to 10*-10° atmos. ; this 
is the order of the instantaneous pressures of 
gaseous products resulting from the detonation 
of solid explosives. Therefore, the fundamental 
importance, from the chemical-physical view- 
point, presented by the knowledge of the proper- 
ties of gases at these particular conditions is 
evident. 
The purpose of this work is in the first place 
to show that, when one admits the validity of 
the thermo-hydrodynamic theory of detonation 





as developed by Chapman,! Jouguet,? and R. 
Becker,’ and makes use of a few experimental 
data,. then one can determine the equation of 
state for gases applicable during the detona- 
tion: the pressures of these gases can be of the 
order of 10° atmos. 

In order to make more comprehensible our 





Fic. 1. Explosion process. 


1—. L. Chapman, Phil. Mag. 47, 90 (1899). 
2 Jouguet, J. de Math. 1, 347 (1905); 2, 6 (1906). 
3 R. Becker, Zeits. f. Physik 8, 321 (1922). 
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procedure, we will sketch briefly the thermo- 
hydrodynamic theory of the explosion. Let us 
consider, for this purpose, Fig. 1, which repre- 
sents, at a certain instant, the state of the 
system composed of the unexploded substance, 
and the products obtained in the detonation, 
when the explosive reaction takes place in a 
cylindrical container with indeformable walls. 
The part of the container on the right J is filled 
with the unexploded substance (gas, liquid, or 
solid), the part on the left JJ is filled by the 
mixture of the explosion products, and the inter- 
mediate part JJI represents the zone (whose 
depth, following Becker, is about 10-7 cm), where 
the explosive transformation takes place. Let us 
denote by D the velocity with which this zone 
(wave front) proceeds (velocity of detonation), 
and W the velocity of the gaseous explosion 
products; let us denote by pi, 71, p1 and by 
p, I, p the density, the temperature, the pressure, 
respectively, of the unexploded substance and of 
the gaseous products. 

In the steady conditions (D=const.), and 
assuming that the effect of heat losses is neg- 
ligible, the principles of continuity of mass, 
conservation of energy, and momentum, then, 
provide the following equations: 


piD=p(D—-W), (1) 


9 


p.D+eD(E.+—) 


(D—W)? 
ness, 6 CR 
2 
pitpD?=pt+p(D—W)? (3) 

(E; and E denote the specific internal energy 
of two phases). 

To the preceding equations, which are hydro- 
dynamical in nature, we shall now add the 
thermodynamic relation : 


= p(D—W)+p(D— w| B+ 


- ¥ 
B-Ey= { cdT-g=0(T-T) 4 (4) 
TI 


and the equation of state for the gases of the 
mixture resulting from explosion : 


F(p, p, T) =0. (S) 


It can easily be shown that Eqs. (1)—(5) permit 
one to obtain D, W, T as functions of the 


parameters ~, p and of the characteristics of the 
unexploded substance : 

















= } 
D=- [P—Pr I. 6) 
Lp—pi pi 
j P—P1 i - 
W=| (p—pi) (7) 
L PPi 
T = ¢(p, p). (8) 
j 
e| 
J 
Vs F K 
ft lite 
Vy vo 


Fic. 2. Hugoniot’s curve for an explosion process. 


Moreover, one arrives at a relation between p 
and v (where v=1/p is the specific volume of the 
gaseous mixture) : 


é.Lo(p, v) —T1]—q=3(p1+p)(%1—2). (9) 


This relation is represented on the (V, p) dia- 
gram (Fig. 2) by a curve known as Hugoniot’s 
curve or H curve. This gives us all the infinite 
states of the gaseous mixture, which are a priori 
possible. In the H curve we can distinguish 
three parts: the first is formed from the complex 
of points above the point G, the second from the 
complex of points between G and F, and the third 
from the complex of points under F. As one can 
immediately verify with the aid of relations (6), 
(7), corresponding to the first part of the H 
curve (detonation) D>W>0; corresponding to 
the second part (which represents no physical 
transformation) D and W assume imaginary 
values; and last, corresponding to the third 
part (deflagration) D>0 and W<0. 

Further, making use of several considerations 
about the entropy of the system, one can prove 
(see for example Schorah‘) that there exists only 
one point of the H curve to which corresponds a 
stable state of the detonation process. This state 
is analytically determined by the fact that, 
corresponding to it, the derivative dS/dv of the 
entropy must be taken equal to zero; graphically 


#R. L. Schorah, J. Chem. Phys. 3, 425 (1935). 
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the stable state is represented by the contact 
point of the tangent to the curve H from the 
point (v1, $1), characteristic of the state of the 
unexploded substance. The fact that, correspond- 
ing to the stable state, we have 


D=U+W 


is of considerable importance, where U denotes 
the velocity of sound in the mixture of gases 
resulting from detonation. 


2 


If we take the equation of state (5) to be 
known, it is possible to calculate all the quanti- 
ties p, p, T, D, W which characterize, from the 
physical standpoint, any gaseous explosive, be- 
cause the pressure of gases deriving from the 
detonation reaches some thousand atmospheres : 
the results of theoretical calculations—see for 
example Lewis'—are in very reasonable agree- 
ment with experimental facts. On the other hand, 
in the detonation of solid explosives gaseous 
mixtures are produced whose pressures can reach 
above 10° atmos. In theoretical calculations one 
experiences the difficulty that at such extremely 
high pressures no experimentally proved equa- 
tion of state is known. In this connection, we 
shall now show that if we start from the experi- 


mental relation between the velocity of detona- 


tion D and the density of the unexploded 


TABLE I. Nitropentaerytrite. 
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substance p:, not only can we unequivocally 
derive the values of the further explosion param- 
eters (W, p, p, T) without a priori knowledge of 
the equation of state for the gaseous products, 
but we can also arrive at such an equation. 
Writing the unknown equation of state in the 
very general form 


pv=r(v)T, (10) 


we see indeed that the H curve of the detonation 
process obeys the equation (f:<)) 


A 


J Seem 


where 


ey 
A=2,T:+g, B(v)=—+H}. 
r(v) 


The stable state, which can be obtained—as we 
have remarked—by the determination of the 
contact point of the tangent to the H curve 
from the initial point (v1, 1), is in this case 
determined by the equation in 7: 


n°’ +2kn—(k—o) =0, (11) 
where 
W 2A (dB/dv) 
n=—, kR=—, o = }n(1—n)vx———(R+7). 
D D? B 


o being very small with respect to k, we can 
easily solve Eq. (11) following a method of 
successive approximations. 

From the values of n, which we obtain in this 






































pl D p p Zz R : 
(kg/dm*)  (m/sec.) (kg /dm‘) (kg/cm?) (°K) way, we derive further the corresponding values 

0.5 3960 0.79 29,600 4670 os 

0.65 4400 1.00 44,900 4840 TaBLE III. Picric acid. 

0.8 4900 1.20 65,000 5050 

1.0 5500 1.45 95,300 5320 - ‘ » r 

t2 6300 1.69 140,500 5720 (kg /dm?) (m /sec.) (kg /dm*) (kg/cm?) (°K) 

1.4 7100 1.93 195,500 6170 

1.6 7900 2.15 262,800 6670 1.03 5150 1.47 83,000 3880 
1.23 5820 1.72 120,700 4070 
1.39 6450 1.93 164,600 4280 

Tapue II. Tetryl. 1.63 7210 2.25 239,400 4660 
D T ; 
(kg /dm®) (m /sec.) (kg /dm*) (kg fem?) (°K) TABLE IV. Trinitrotoluene. 

1.00 5480 1.43 91,800. 4400 » . a 

1.28 6510 1.80 160,400 4740 o P ‘ 

1.45 7220 2.02 218,100 4980 SS ES a dha 

1.54 7375 2.15 242,500 5100 1.0 4700 1.44 68.700 3210 

1.61 7470 2.20 259,100 5140 1.29 5900 1.82 132,800 3610 
1.46 6500 2.04 178,000 3860 
1.59 6900 2.22 216,200 4020 


5 Lewis, Chem. Rev. 10, 49 (1932). 
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Fic. 3. Equation of state for detonation gases: 
nitropentaery trite. 
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Fic. 4. Equation of state for detonation gases: tetryl. 


of p, p, T making use of Eqs. (6)—(8). Following 
this procedure we have the possibility of establish- 
ing the function r(v) and therefore the unknown 
equation of state. 

We have performed numerical computations 
for several explosives (nitropentaerytrite, tetryl, 
picric acid, trinitrotoluene) assuming as starting 
point the experimental relation between the 
detonation velocity D and the density of explo- 
sive, derived from the results of careful measure- 
ments performed by Friederich.* For the chemical 
equations of the explosive reactions, character- 
istic of the detonation processes, we have taken 
the data from A. Schmidt’; the values of the 
mean specific heats as functions of the tempera- 
ture, for every gaseous mixture were calculated 
starting from the data contained in Beyling- 
Drekopf’s book.® 

Tables I-IV present a summary of the results 
of our numerical calculations. 

From our results we can conclude that the 
quantities p,v, T of the gaseous mixtures re- 


6 W. Friederich, Zeits. ges. Schiess. und Sprengstoffe 28, 
2, 51, 80, 143, 213, 244 (1933). 

7A. Schmidt, Zeits. ges. Schiess. und Sprengstoffe 33, 
364 (1938); 34, 8, 37, 80 (1939). 

§ C. Beyling and K. Drekopf, Sprengstoffe und Ziindmittel 
(Julius Springer, Berlin, 1936). 
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Fic. 5. Equation of state for detonation gases: picric acid. 
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Fic. 6. Equation of state for detonation gases: 
trinitrotoluene. 


sulting from the several processes of detonation 
are in good accordance with the equation (due to 





Van Laar) 
plv—a(v) ]=nRT, (12) 
with 5 
a(v) = iat ’ (13) 
1+c/v 1+¢p 


where the constants b,c have the following nu- 
merical values: 


nitropentaerytrite: b=2.08 c=2.00 
tetryl: b=1.81 c=1.62 
picric acid: b=1.96 c=1.79 
trinitrotoluene: b=2.07 c=1.88 


In order to illustrate the accuracy with which the 
equation of state (12) is verified, we have plotted 
in Figs. 3-6 the values of p against the values of 
1/a for the gaseous mixture of the explosives 
considered. One can see that the dependence is 
really linear, as one must expect from Eq. (13). 

Finally, attention is called to the fact that the 
equation of state (12), which we. have derived, 
is suitable only within the range of pressures and 
temperatures occurring in processes of detonation 
of solid explosives, that is approximatively from 
30,000 to 300,000 atmos., and from 3000 to 
7000°K. 
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The Critical Micelle Concentration 
of Anionic Soap Mixtures 


H. B. KLEVENS* 


Department of Chemistry, University of Chicago, 
hicago, Illinois 


September 5, 1946 


ANY technical papers have dealt with the increase in 
solubilization, wetting, detergency, etc., of soap 
mixtures as compared with pure soaps, but except for some 
ultramicroscopic studies there have been no data reported 
on the properties of soap mixtures in the region of the 
critical micelle concentration (cmc). The changes occurring 
in dye spectra upon the addition of varying amounts of 
soaps! have been used in this work. 

The change in cmc (the sum of the concentration of the 
two soaps) of soap mixtures (potassium myristate plus 
potassium laurate, potassium caprate, potassium capry- 
late, and sodium lauryl sulfate) as determined with the 
dye, pinacyanol chloride, is plotted in Fig. 1. The cmc of 
potassium caprylate and pdtassium caprate have not been 
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Fic. 1. The critical micelle concentrations in moles per liter (the 
sum of the concentrations at the appearance of micelles) of various 
soap mixtures. 
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reported previously by means of the dye method, but the 
values of 3.94 10~! and 9.9 10-? M observed here agree 
with those of 3.6—3.8 X 107 and 9.5—10.0 x 10- reported by 
Ekwall,? and 3.7-4.0X10-! and 10.0-10.210-? as de- 
termined by refraction.’ Other unreported data which will 
be discussed in a following report include mixtures of 
alkyl sulfonates as well as potassium laurate plus potassium 
caprate and potassium caprylate, and potassium caprate 
plus potassium caprylate. Similar results have been ob- 
served for mixtures of cationic soaps using indophenol‘ as 
the dye indicator. 

It can be seen from the curves in Fig. 1 that the greatest 
change in cmc of the soap mixtures is observed when there 
is the largest difference in the cmc of the pure soaps. It is 
also evident that there are no cmc values which do not fall 
within the range of values of the pure soaps. These results 
over certain mole fraction ranges are similar to the effects 
of salts on the cmc of pure soaps.® In this connection, how- 
ever, it is necessary to consider not the sum of the concen- 
tration in the soap mixture at the appearance of micelles, 
but the concentration of one of the soaps at this point as 
affected by the other. When this is done, it is found, for 
example, that potassium caprate acts like a salt, say 
potassium chloride or sodium chloride, depressing the 
necessary concentration of potassium laurate. This value, 
which is equal to the mole fraction times the total cmc, is 
found to be similar to a salt effect for mole fractions of the 
laurate as low as 0.50. 

This treatment was found to be valid in various other 
mixtures only in the region of higher mole fractions of the 
less soluble soap. In the portion of the curves in Fig. 1 
where the greatest change in the total cmc occurs, the salt 
effect phenomena were not observed. Thus, a much greater 
change in the total cmc can be obtained by adding a small 
amount of a less soluble soap than by adding an equivalent 
amount of a salt. The tendency of the less soluble soap to 
aggregate seems to have a much greater effect on the total 
aggregation associated with micelle formation than a salt. 

When two soaps of approximately the same cmc are 
mixed, there is no observable change in the cmc of the 
mixture and no region of a salt effect is observed. A mixture 
of this type, potassium myristate plus sodium lauryl 
sulfate, is seen in Fig. 1. 

In the consideration of soap mixtures, it seems that the 
difference in the tendency of single soaps to form micelles 
controls the micelle formation in their mixtures; the greater 
the difference (for example in solubility), the greater will 
be the initial effect of the more soluble soap. If this tend- 
ency to aggregate is different in two soaps, then in a 
limited concentration range of mixtures, the soap having 
the lower tendency to aggregate will act as a salt towards 
the other. Soaps having about the same tendency to aggre- 
gate will have little effect on each other as far as salt effect 
and change in total cmc are concerned. 

* Present address: Chemical & Physical Research Laboratories Fire- 
stone Tire & Rubber Company, Akron, Ohio. 

1M. L. Corrin, H. B. Klevens, and W. D. Harkins, J. Chem. Phys. 
14, 480 (1946). 

2P. Ekwall, Kolloid Zeits. 101, 135 (1942). 

3H. B. Klevens, J. Chem. Phys. 14, 567 (1946). 


4H. B. Klevens, J. Am. Chem. Soc. to be published. 
5M. L. Corrin and W. D. Harkins, in preparation. 
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A Relation Between Bond Multiplicity and 
Interatomic Distance 


R. T. LAGEMANN 
Emory University, Emory University, Georgia 
November 14, 1946 


N a note bearing the above title, Kavanau! has ad- 
vanced a most useful relation connecting the interatomic 
distances of covalent bonds to the bond multiplicity. He 
suggests that 
aan) (1) 
m+n2+1 


where D is the interatomic distance in angstroms; @ and b 
are constants determined by the specific atom pair; N is 
the bond multiplicity; and m, and m2 are the principal 
quantum numbers of the valence electrons in the bonded 
atoms. 

Examination of the constants a and b given by Kavanau 
shows them to be roughly related to the sum Z of the 
atomic numbers of the atoms making up the bonded pair. 
If we set b proportional to Z and a inversely proportional 
to Z, we have 


D=a+0( 











Ki ny +n —1 N 
=F +Ke(“), (2) 
PA m+n2+1 
TABLE 1. Observed and calculated bond lengths. 
Bond Bond 
length length 

Atom Bond in A in A Percent 

Molecule pair order D obs. Decale. deviation 
C2He cc 1.11 1.55¢ 1.50 3.2 
CH;CCH CC Lae 1.46¢ 1.47 0.7 
CeHe cc 1.78 1.39¢ 1.36 pe 
4 cx 2.12 1.33¢ 1.31 1.5 
CH;:CCH cx 2.89 1.21¢ 1.22 0.8 
C2He CC 3.00 1.2048 1.21 0.8 
CH:3F CF 1.00 1.3969 1.44 2.9 
CHsNH2 CN 1.00 1.478 1.49 1.4 
N CN 3.00 1.1544 1.15 0.0 
CH;0H co 1.00 1.438 1.46 21 
H:CO co 2.00 1,219 1.23 iP 
Fe FF 1.00 1.4354 1.42 1.4 
NH2N He NN 1.00 1,.47> 1.46 0.7 
(CeHsN )2 NN 2.00 1.236 1,23 0.0 
Ne NN 3.00 1.095* 1.09 1.0 
OF 2 OF 1.00 1.418 1.42 0.7 


_ 
w 


Average percent deviation 








a W. Gordy, J. Chem. Phys. 14, 305 (1946). 

b J. L. Kavanau, J. Chem. Phys. 12, 467 (1944). 

eR. S. Mulliken, C. A. Rieke, and W. G. Brown, J. Am. Chem. Soc. 
63, 41 (1941). 


where K; and K2 are considered constant for all pairs of 
atoms having the same value of m and the same value of 
nz. This modification covers more cases with two constants 
than does Kavanau’s original relation, and seems to have 
considerable merit in empirically correlating bond length 
and multiplicity. Figure 1 shows the linearity of DZ versus 
Z* (3/5)* for pairs of atoms for which m =n2=2. 

Table I gives a list of atom pairs whose valence electrons 





24 


Oz 
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possess principal quantum numbers of 2. Included in the 
table are values of the observed interatomic distances and 
values of the interatomic distance calculated from Eq. (2) 
using K,=12.45, and K:=0.0676. The values of K, and 
Kz were obtained by a least squares treatment of the ex- 
perimental bond lengths. The combination m,=n2.=2 was 
selected at random for presentation here, but other com- 
binations of the principal quantum numbers appear on 
preliminary inspection to yield similar results. 

It will be noted that the proposed relation holds for 
isotopic molecules where Z and D are constant. However, 
it does not care in a precise way for cases such as CO and 
Ne where Z is the same for both molecules but the observed 
interatomic distance is clearly not. 


1J. L. Kavanau, J. Chem. Phys, 12, 467 (1944). 
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